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Abstract 


Both space and laboratory plasmas can be associated with static magnetic field, and the 
field geometry varies from uniform to non-uniform. This thesis investigates the impact 
of magnetic geometry on wave modes in cylindrical plasmas. The cylindrical configura¬ 
tion is chosen so as to explore this impact in a tractable but experimentally realisable 
configuration. Three magnetic geometries are considered: uniform, focused and rippled. 

For a uniform magnetic field, wave oscillations in a plasma cylinder with axial flow and 
azimuthal rotation are modelled through a two-fluid flowing plasma model. The model 
provides a qualitatively consistent description of the plasma configuration on a Radio 
Frequency (RF) generated linear magnetised plasma (WOMBAT, Waves On Magnetised 
Beams And Turbulence [Boswell and Porteous, Appl. Phys. Lett. 50, 1130 (1987)]), and 
yields agreement between measured and predicted dependences of the wave oscillation 
frequency with axial field strength. The radial profile of the density perturbation predicted 
by this model is consistent with the data. Parameter scans show that the dispersion curve 
is sensitive to the axial field strength and the electron temperature, and the dependence of 
the oscillation frequency with electron temperature matches the experiment. These results 
consolidate earlier claims that the density and floating potential oscillations are a resistive 
drift mode, driven by the density gradient. This, to our knowledge, is the first detailed 
physics modelling of plasma flows in the diffusion region away from the RF source. 

For a focused magnetic field, wave propagations in a pinched plasma (MAGPIE, MAG- 
netised Plasma Interaction Experiment [Blackwell et ah, Plasma Sources Sci. Technol. 
21, 055033 (2012)]) are modelled through an ElectroMagnetic Solver (EMS) based on 
Maxwell’s equations and a cold plasma dielectric tensor.[Chen et. ah, Phys. Plasmas 
13, 123507 (2006)] The solver produces axial and radial profiles of wave magnitude and 
phase that are consistent with measurements, for an enhancement factor of 9.5 to the 
electron-ion Coulomb collision frequency and a 12% reduction in the antenna radius. It 
is found that helicon waves have weaker attenuation away from the antenna in a focused 
field compared to a uniform field. This may be consistent with observations of increased 
ionisation efficiency and plasma production in a non-uniform field. The relationship be¬ 
tween plasma density, static magnetic field strength and axial wavelength agrees well with 
a simple theory developed previously. Moreover, the wave amplitude is lowered and the 
power deposited into the core plasma decreases as the enhancement factor to the electron- 
ion Coulomb collision frequency increases, possibly due to the stronger edge heating for 
higher collision frequencies. 

For a rippled magnetic field, the spectra of radially localised helicon (RLH) waves 
[Breizman and Arefiev, Phys. Rev. Lett. 84, 3863 (2000)] and shear Alfven waves (SAW) 
in a cold plasma cylinder are investigated. A gap-mode analysis of the RLH waves is 



Xll 


first derived and then generalised to ion cyclotron range of frequencies for SAW. The 
EMS is employed to model the spectral gap and gap eigenmode. For both the RLH 
waves and SAW, it is demonstrated that the computed gap frequency and gap width 
agree well with the theoretical analysis, and a discrete eigenmode is formed inside the 
gap by introducing a defect to the system’s periodicity. The axial wavelength of the gap 
eigenmode is close to twice the system’s periodicity, which is consistent with Bragg’s law, 
and the decay length agrees well with the analytical estimate. Experimental realisation 
of a gap eigenmode on a linear plasma device such as the LArge Plasma Device (LAPD) 
[Gekelman et ah, Rev. Sci. Instrum. 62, 2875 (1991)] may be possible by introducing 
a symmetry-breaking defect to the system’s periodicity. Such basic science studies could 
provide the possibility to accelerate the science of gap mode formation and mode drive in 
toroidal fusion plasmas, where gap modes are introduced by symmetry-breaking due to 
toroidicity, plasma ellipticity and higher order shaping effects. 

These studies suggest suppressing drift waves in a uniformly magnetised plasma by 
increasing the field strength, enhancing the efficiency of helicon wave production of plasma 
by using a focused magnetic field, and forming a gap eigenmode on a linear plasma device 
by introducing a local defect to the system’s periodicity, which is useful for understanding 
the gap-mode formation and interaction with energetic particles in fusion plasmas. 
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predictions for the gap centre (xo) and the lower and upper edges of the gap 
(x_ and X-), (b) spatial distributions of wave energy for antenna frequencies 
inside and outside the spectral gap (solid bar denotes the antenna location). [83] 
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Chapter 1 


Introduction 


“A plasma is a quasineutral gas of charged and neutral particles which exhibits collec¬ 
tive behaviour” [12] It is often referred to as the “fourth” state of matter, following the 
well-known solid, liquid and gaseous states, because of the special collective behaviour 
involved. It was first identified by Sir William Crookes in 1879 in a Crookes tube, who 
named it “radiant matter”. m The term “plasma” was coined by Irving Langmuir in 1927 
considering that the electrified fluid carries high velocity electrons, ions and impurities in a 
similar way of the blood plasma which carries red and white corpuscles and germs. mm 
The Greek word “7rAd(x//a” (“plasma”) means “moldable substance”. Since the mercury 
arc that Irving Langmuir used to study oscillations in ionised gases diffuses through the 
whole glass chamber, and molds itself, this may be another reason for Irving Langmuir to 
call the ionised gas “plasma”. m 

Although the plasma rarely exists in our natural life due to its high temperature, more 
than 90% of the matter of the universe is believed to be composed of plasmas, e. g. stellar 
interiors and atmospheres, gaseous nebulae and much of the interstellar medium, excluding 
the more speculative nature of dark matter. nacET] The field of plasma physics can be dated 
back to 1920s when Langmuir, Tonks and their collaborators worked on gas discharges. It 
was significantly advanced later by the nuclear fusion program and industrial applications 
of plasma processing. mm Plasma physics is the underlying science of space and solar 
physics, astrophysics, and finds applications in diverse areas of MagnetoHydroDynamics 
(MHD) energy conversion and ion propulsion, solid state plasmas and gas lasers. [12] 

Both space and laboratory plasmas can be associated with static magnetic field, and 
the field geometry varies from uniform to non-uniform. This thesis investigates the impact 
of magnetic geometry on wave modes in cylindrical plasmas. The cylindrical configura¬ 
tion is chosen so as to explore this impact in a tractable but experimentally realisable 
configuration. Three magnetic geometries are considered: uniform, focused and rippled. 

This chapter first introduces the plasma source employed in the thesis, and then il¬ 
lustrates four magnetic geometries. Specifically, the dispersion relation of helicon waves, 
which drive the helicon discharge and will be considered in detail in the following chapters, 
is derived from a cold-plasma dielectric tensor perspective. Then the concept of helicon 
discharge is given, together with its typical applications. We select low temperature helicon 
plasma sources for comparison to experiments, both because of availability, and diagnostic 
opportunities. Finally, examples of helicon discharge with uniform (WOMBAT, Waves On 
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Magnetised Beams And Turbulence[20j) and focused magnetic fields (MAGPIE, MAGne- 
tised Plasma Interaction Experiment [8]), large linear plasma device with multiple magnetic 
mirrors (LAPD, LArge Plasma Device^]), and magnetic confinement fusion with toroidal 
magnetic field (tokamak) are given. Gap eigenmodes that will be formed in a linear plasma 
with rippled magnetic held, by introducing a local defect to the system’s periodicity, are 
applicable to understanding the gap-mode formation and its interaction with energetic 
particles in tokamaks. ]2T] 


1.1 Plasma source 

Although a cold plasma model, which assumes zero-temperature frictionless fluids of ions 
and electrons, ignores physics such as pressure and how, it is tractable in its mathematical 
analysis and provides a remarkably accurate description of the common small-amplitude 
perturbations that are possible for a hot plasma. [22] Particularly, it is applicable to plasmas 
and physics for which the thermal effects and particle motions are ignorable, e. g low- 
temperature plasma and long-time scale problems. For the helicon discharge studied in 
this thesis, the temperature is low, and so we utilise the cold plasma model to describe 
helicon waves and shear Alfven waves. The approach that is followed is essentially that of 
Stix, and Gurnett and Bhattacharjee. [22l [23] . 


1.1.1 Helicon wave 

Conductivity tensor and dielectric tensor 

We start from the hrst of Maxwell’s equations 

„ „ . dE dD , . 

VxB = /j 0 j + /ioeo-^- = /io^-, (1.1) 

where the plasma current j and the vacuum displacement eo<9E /dt have been combined 
into the electric displacement D. Here, E and B are the electric and magnetic fields, 
respectively, //o is the permeability of free space, eo is the permittivity of free space, and 
t is time. The conductivity tensor d and the dielectric tensor e are defined as: 


j = o- ■ E, D = e 0 e E, 


( 1 . 2 ) 


respectively. We substitute Eq. (1.2) into Eq. (1.1) and consider perturbations of the 
general form exp[*(k • r — cut)], with k the wave vector, r the space vector, and u ; the wave 
frequency. The relation between a and e is thus 


e = I-—, (1.3) 

iuje o 

with 1 the unit tensor and i the unit imaginary number. The plasma is assumed to be 
uniform, unbounded and immersed in a uniform static magnetic field Bq. For a cold 
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plasma, the equilibrium velocities of particles are zero, and the equilibrium electric field 
vanishes. Therefore, a set of linear equations that govern the plasma system is formed: 


E = Ei, B = B 0 + Bi, v s = v s i, n s = n s0 + n s i, 


(1.4) 


where v s is the particle velocity and n s is the plasma density. The subscript s denotes 
the plasma species, namely electrons and ions, 0 denotes the equilibrium values, and 1 
denotes the small-amplitude perturbations. The linearised current density and momentum 
equation are thus: 

jl = X^s^sOVsl, (1-5) 


3v s i „ 

= q s El + q s v s 1 x B 0 , 
ot 


( 1 . 6 ) 


with q s and m s the particle charge and mass respectively. Substituting the solved v s i 


from Eq. (1.6) into Eq. (1.5) and recalling ji = a ■ Ei, we get the conductivity tensor 


<7 = 


n s0 q s 


m. 


bJcs 


—iuj 




0 

0 


0 


0 


(1.7) 


where uj cs = q s \Ro\/m s is the cyclotron frequency. According to Eq. (1.3), the dielectric 
tensor is thereby [22] 

S -iD 0 " 

e = %D S 0 . (1.8) 

OOP 

Here, the quantities S (for sum), D (for difference), and P (for plasma) are defined as: 


S = (R + L)/ 2, D = (R-L)/ 2, 


(1.9) 


R = 1-E 


uj ; 


ps 


l = 1 -E 


CO, 


ps 


> ’ ^ a; 2 ’ 


( 1 . 10 ) 


UJ(u+U) cs ) , s w(w - Ucs) : 

and u)p S = \/q%n s o/{eQm s ) is the plasma frequency. It will be shown later (through 
Eq. (1.17) and Eq. ( 1.18| >) that the parameters R and L stand for right-hand circularly 
polarised (RHCP) and left-hand circularly polarised (LHCP) modes, respectively. Equa¬ 


tions (1.8)- (1.10) represent the well-known cold-plasma dielectric tensor which has been 


used widely for plasma wave analyses, and will be employed below to derive the dispersion 
relation of helicon waves. 


Whistlers: unbounded helicon waves 


We first study the normal modes of this unbounded plasma system. Combining Eq. (1.1) 
with the Faraday’s law 


V x Ei = - 


0B X 

~df' 


( 1 . 11 ) 
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we get the homogeneous-plasma wave equation 


co 


k x (k x Ei) + —^e • Ei = 0, 


( 1 . 12 ) 


which can be written in a more convenient form 


n x (n x Ei) + e • Ei = 0 


(1.13) 


by introducing the dimensionless index of refraction vector n = k c/co. The magnitude of 
n = |n| = c/ (co/k) is the ratio of the speed of light c = 1/ y/Jipeo to the wave phase velocity 
Uph = oo/k, with k = |k| the wave number. Without losing the generality, we choose a 
3D Cartesian coordinate system ( x , y , z ) and assume n to be in the x — z plane and Bq 


along the z axis. The angle between n and Bo is d. Equation (1.13) can be written in a 
matrix form [22] 


S — n 2 cos 2 id —iD n 2 sin id cos d 
iD S-n 2 0 

n 2 sin d cos d 0 P — n 2 sin 2 d 



E\± 


E\y 


E\z 


= 0 , 


( 1 . 14 ) 


which has nontrivial solutions of Ei only if the determinant of the matrix equals zero. 
The zero-value determinant results in 


tan 2 d ~ P (n 2 ~ R)(n 2 - L) 
(Sn 2 - RL){n 2 - P) ’ 


(1.15) 


describing the normal modes (no source term in Eq. (1.13)) of the system. For modes prop¬ 


agating parallel to the static magnetic field, namely d = 0, the electric field eigenvectors 


associated with each of the three roots of Eq. (1.15) are: 


l — ( 0 , 0 , Eiz), 

( 1 . 16 ) 

= (Eixi 'i / E\x') 0 ), 

( 1 . 17 ) 

“ (Eixi 0 ). 

( 1 . 18 ) 


Here, the y component of the electric field has been expressed in form of the x component. 


Equation (1.16) represents the electron plasma oscillations which occur at the plasma 
frequency according to Eq. (1.10), i. e. co 2 = w 2 j + co 2 


pe 


co 2 e ( m e <C mi), and do not 


propagate. This is the same as that in an unmagnetised plasma because the presence of 
the static magnetic field does not affect particle motions (driven by ) in the parallel 
direction. This mode is called longitudinal mode due to the parallel wave vector to the 


electric field. Equations (1.17) and (1.18) represent transverse modes because of k T Ei, 


and from the phases of E\x and E\y we can see that they are RHCP and LHCP modes, 
respectively. The dispersion relations of these two modes are plotted in Fig. |l.l|for a singly 


ionised argon plasma, as an example, with density n s o = 10 18 m 3 and field strength 
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|Bo | = 0.01 T. While Fig. 0 is the straight visualisation of Eqs. (|1.17|)-(p~.18|), as done 



Figure 1.1: Dispersion curves of waves propagating parallel to the static magnetic field in a 
cold, unbounded, and uniform plasma: I, shear (or slow) Alfven waves; II, compressional (or fast) 
Alfven waves which is also called whistler mode; III and IV are free space modes for w > ojl (cut-off 
frequency for the LHCP mode) and u> > uir (cut-off frequency for the RHCP mode), respectively. 


by Ginzburg, |21J the singularity at uj = 0 is “unreal” and can be removed by introducing 
the quasineutral condition of q e n e o + qiUio = 0. In the fact, for this low frequency range 
(uj <C Uci) of waves, careful analysis shows that[23] 

n 2 « 1 -I— 2 * = 1 + —o = n A- (1.19) 

“ci V A 

Here, va = Bq / y/fJo r Ho m i and ua are the phase velocity and refraction of index of Alfven 
waves, as labeled in Fig. 0 Branch I is called shear (or slow) Alfven waves (SAW), which 
will be used to form a gap eigenmode in Chapter [5] Due to the same direction of the ion 
cyclotron motion and wave polarisation, this branch sees the ion cyclotron resonance 
when uj -A- u c i , forming the so-called ion cyclotron range of frequency (ICRF) waves. 
The ICRF waves have an important role in heating ions in fusion plasmas. [25] Similarly, 
the branch II, named compressional (or fast) Alfven waves, sees the electron cyclotron 
resonance when uj -A- |w ce |, because of the same direction of the electron cyclotron motion 
and wave polarisation. The electron cyclotron range of frequency (ECRF) waves have 
been used widely in studying the energy transport, driving the current, and diagnosing 
the radial profile of electron temperature. Waves polarised in the opposite direction of 
particle cyclotron motions do not resonate with these particles, so that the LHCP branch 
does not see the electron cyclotron motion and the RHCP branch does not see the ion 
cyclotron motion. Branch II is also called whistlers, which for uj c i uj <C \uj C e\ are viewed 
as unbounded helicon waves. m The dispersion relation of whistlers for uj c i <C uj <C \uj ce \ 
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can be obtained by neglecting the ion terms and assuming \oo ce \ <C oo pe in Eq. (1.17) 

,2 


n 2 = 


oo, 


pe 


UJ, 


pe 


u(\u ce \-u) w|a lce\' 


( 1 . 20 ) 


which further gives the corresponding group velocity 


oo 


Vg = 2C- 


1/2 (I 


00 r . 


~^) 3/2 -, 2 ^ 1 / 2 I ^ I 1/2 


00 . 


( 1 . 21 ) 


pe 


Equation (1.21) shows that the group velocity of whistlers increases with the wave fre¬ 


quency. This implies that a received frequency of a finite duration burst will decrease in 
tone, leading to the name “whistler”. It was first observed around the second half of the 
World War I,[IS] reported later by Barkhausen in 1919. [26l [27] and explained by Storey 
in 1953.|28j It can be used to measure the electron density due to n 2 oc n e , and believed 
to be the only diagnostic for the electron density at high altitudes above the Earth’s iono¬ 
sphere, before direct in situ measurements were available. [23J Interestingly, there exists 
another type of whistlers which, however, have a rising tone with the increasing time. 
They occur near the ion cyclotron frequency on the branch I, and are thereby called ion 
cyclotron whistlers. [29] The ion cyclotron whistlers cannot be observed from the ground 
but satellites outside of the ionosphere. [ 


Helicon waves in a cylindrical uniform plasma 


When the plasma is bounded, e. g. in a cylindrical geometry, the wave field of whistlers 
needs to meet the either conducting or insulating boundary conditions, and their pure 
electromagnetic features cannot be kept. Then they become helicon waves, named after the 
associated helical force lines which rotate and carry electrons.] 18] The word “helicon” was 
first suggested by Aigrain in 1960 to describe an electromagnetic wave which propagates in 
solid metals at low temperatures and with frequencies oo C i <C oo <C |w ce |. m The dispersion 
relation of helicon waves in a cylindrical plasma with uniform density was first derived by 
KMT, [31] and more recently by Chen and Arnush. [5211551151] We shall follow the steps of 
the latter. 

For a cylindrical coordinate system (r, 6, z), the exponential factor changes to 
exp [i(m0 + k z z — oot)] with k z the axial wave number and m the azimuthal mode number. 
We introduce a total current jt = dT)/dt and recall D = e^e • E to get 



s 

RL 

iD 

RL 

0 ' 


jtlr 

hje 0 Ei = e 1 • jti = 

iD 

RL 

S 

RL 

0 


jtie 


0 

0 

1 

p . 


jtlz 


( 1 . 22 ) 


which can be equivalently written as 


iwe 0 Ei = ajti + ia h ( z x j t i) + a d [i x (z x jti)]- 


(1.23) 
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Here, the relation S 2 — D 2 = RL (according to Eq. (1.9)) has been used, a c = 1/P is 
the conduction or polarisation current, ah = —D/(RL) is the Hall current or E x B 
drift, ad = 1/P — S/(RL ) is the displacement current, and z is the unit vector along 
z.|33j Assuming that the plasma is sufficiently dense (|o; ce | <C co pe ) to ensure the ignorable 
displacement current compared to the plasma current and the wave frequency is in range of 


i/w c i|w ce | « a; « |w ce | to ignore the ion current, we neglect the ion terms in Eqs. (ILlOj) 
and get 


, ,2 , ,2 

U/np UJrtp 

__ P e T _ P e 


+ W ce ) 


, L = - 


u(u - u ce ) 


p _ u, p e 
OJ 2 • 


(1.24) 


Equation (1.24) results in a c = — ui 2 /uj 2 ah = -wwce/wi, ad = 0, and simplifies 


Eq. (1.23) to 


Ei = 


Ur. 


e ou pe 


(iSrtj tl + Z X j t i) 


(1.25) 


with 5 rt = —u/u.) ce . Combining Eq. (1.25) with Eq. 0 (neglect the vacuum displacement 
eo<9Ei /dt) and Eq. (1.11), we have 

5 r tS7 x V x Bi — k z V x Bi + fc^Bi = 0 (1.26) 

with k/u = — poujn e oq e /B q = 5 r tk 2 . Here, k w is the wave number of helicon waves in 


free space (namely whistlers according to Eq. (1.20)), and k s = u pe /c is the skin number 


standing for the decay constant of electromagnetic waves penetrating into the plasma. 


Equation (1.26) can be factored into[31 j 


(A-Vx)(ft-Vx)B 1 = 0, 

with /3i and /?2 the separation constants representing the roots of 

brtfi 2 — k z /3 + k^, = 0. 


(1.27) 


(1.28) 


In the limit of m e = 0, which yields S r t = 0, there is only one root of Eq. (1.28), 

(d fJ'O^'eOQe 


k 2 


kz Ho 


(1.29) 


This is the usual dispersion relation of helicon waves in a bounded uniform plasma. The 
root /3, as will be seen in Eq. (1.36), is determined by the boundary condition and the 


mode structure of Bi. Taking the curl of Eq. (1.26) and employing V • Bi = 0, we get 

V 2 Bi +/3 2 Bi = 0. (1.30) 


The axial component of Eq. (1.30) is 


d 2 B lz | 1 dB lz | 




B lz = 0 , 


(1.31) 
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which has the solution 


Biz = AJ m (Tr) 


(1.32) 


due to the finite B\ z at the origin. Here, A is the amplitude constant, J m is the Bessel 
function of the first kind of order m, and T 2 = (3 2 — k 2 . The other two components of Bi 


can be obtained from Eq. (1.26) (5 r t = 0): 

+ k z )J m -i(Tr) + (/3- k z )J m+ i(Tr)\, 


Bir — 


iA , 


2 T L 
A 


Big = - — [(/3 + k z )J m -i(Tr) -((3- k z )J m+ i(Tr) 


(1.33) 


(1.34) 


Equations (1.32)-(1.34) describe the radial and azimuthal components of the wave mag¬ 


netic field, from which we can also compute the electric wave field through Eq. (1.11): 


E\r — r—BiQ, E\o — — — Bir, Ei z — 0. 
k z k z 


(1.35) 


Here, the axial component of E vanishes because of m e = 0 assumed in Eq. (1.25), namely 
E\ z = —{iujm e /q 2 n e o)jtiz = 0. A conducting boundary requires Eig(a) = Ei z (a) = 0 
which thereby gives Bi r (a) = 0, while an insulating boundary requires ji r (a) = 0. Here, a 
is the radius of the boundary wall. Neglect the displacement current in Eq. (1.1) (w/c 2 -C 


1) and combine it with Eq. (1.26) (5 r t = 0), we get ji = (/3//ro)Bi which implies a 


vanishing radial component of B, i. e. Bi r (a) = (no/fi)jir{p) = 0. Hence, there exists a 
universal boundary condition for both conducting and insulating boundaries for m e = 0, 
namely, 

(/3 + k z )J m -i(Ta) + (/3- k z )J m+ i(Ta) = 0. (1.36) 


The root /3 can then be solved from Eq. (1.36) for a given mode structure. For the 
lowest two azimuthal modes, i. e. m = 0 and m = 1 which are usually observed in 


helicon discharges, Eq. (1.36) gives Ji(Ta) = 0 and Ji(Ta) = (Tak z /2f3)[J2(Ta) — Jo (To)], 
respectively, which also vanishes for ak z C 1, a long, thin tube. The lowest root of 
Ji(Ta) = 0 is Ta = 3.83 which yields /3 ~ 3.83/a for k z <S 1/a. To illustrate the helicon 
wave, we shall plot the electromagnetic field line patterns of the m = 0 and m = 1 modes. 
Putting back the exponential factor, we can write the real parts of the wave field of these 

two modes as: 

(m = 0 ) 

B\ r (r,z,t) = j;k z Ji(XV) sin(k z z — u>t), 

Big(r, z, t ) = ji/3Ji(Tr) cos(k z z — cot), 

Bi z (r, z, t ) = AJo(Tr ) cos(k~z — cot), 

Ei r (r,z,t) = y ^/3Ji(Tr) cos(k z z - ut), 

Eig(r, z, t ) = —^uJi(Tr) sin (k z z — cut), 

Ei z (r,z,t) = 0, 


(1.37) 
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(m = 1 ) 

B\r = —£r[(P + k z )J 0 (Tr) + (/3 - k z )J 2 (Tr)\ sin(6> + k z z - cot), 

Bio = + k z )Jo(Tr) - (/3 - k z )J 2 {Tr )] cos(0 + k z z - Lot), 

B\ z = AJ\(Tr) cos(9 + k z z — ot), (1.38) 

Eir = ~4rTt\(.P + k z )J 0 (Tr) - (/3 - k z )J 2 (Tr )] cos(0 + k z z - Lot), 

Eio = ^tjr[(P + k z )J 0 (Tr) + (/? - k z )J 2 {Tr)\ sin (6 + k z z - ut), 

Eiz = 0. 


The cross sections of these two modes can be visualised easily by comparing the radial 


and azimuthal components in Eq. (1.37) and Eq. (1.38). For m = 0, we introduce @o = 
B\ r /Big, an d plot the electromagnetic field line pattern as a function of qq. As shown in 
Fig. |1.2[ the magnetic field line pattern changes from purely radial to spiral and to purely 
circular as |go| is decreased. The electric field line pattern shows a reverse trend, consistent 



Figure 1.2: Electromagnetic field line patterns for the m = 0 mode as a function of g 0 where 
go = Bir/Bio: upper row shows the magnetic held line pattern and lower row shows the electric 
field line pattern. The axes represent the cross section edges of the plasma cylinder. The patterns 
for go = 0 represent a helicon mode. 

with Faraday’s law. Due to the usually true \k z /{3\ <C 1 in helicon discharges (inclusion of 
cylindrical boundary decreases k z over the case of infinite plane wave propagation), [36j 32] 
which implies | go I 1 an d the dominance of E\ r over E\q, the best way to excite this mode 
is through coupling the E\ r component. Here, E\ r and E\q stand for the electrostatic and 
electromagnetic components, respectively. For m = 1, we introduce q\ = (1 + k z /j 3)/( 1 — 
k z /(3) and set k z z — uot = 0. Figure [L3] shows the corresponding magnetic and electric field 
line patterns which stay almost the same with g\ increased from 1/3 to 9/11, 1, 11/9 and 
3. The corresponding values of k z //3 are: —1/2, —1/10, 0, 1/10 and 1/2, respectively. As 
time advances or z increases, the field line pattern rotates about its centre point but the 
shape does not change. To show the variation of wave field strength in radius, we plot the 
wave magnetic field over radius for the m = 0 and m = 1 modes. Here, all components are 
normalised to their own maximum values to show a clear radial structure. Figure |l~4| ]a.) 
shows the result of the m = 0 mode, which implies that the profiles of B\ r and B\g are 
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= 1/3 0 i = 9/11 0 i = 1 0 i = H/9 0 i = 3 



Figure 1.3: Electromagnetic field line patterns for the m = 1 mode as a function of qi where 
0 i = (1 + k z //3)/( 1 — k z //3): upper row shows the magnetic field line pattern and lower row shows 
the electric field line pattern. The axes represent the cross section edges of the plasma cylinder. 
The patterns for 0 i = 1 represent a helicon mode. 


overlapped and peak off axis, while the profile of B\ z peaks on axis. By contrast, for the 


(a) m = 0 



( b ) m = 1 



Figure 1.4: Radial structures of the m = 0 mode (a) and m = 1 mode (b). 


m = 1 mode as shown in Fig. |1.4[ b), the profiles of B\ r and B\g peak on axis, while the 
profile of B\ z peaks off axis and vanishes at the origin. The difference between Fig. |1.4[ a) 
and Fig. |1.4[b) helps identifying m = 0 and m = 1 modes in experiments. 


Although the approximation of m e = 0 works for most scenarios of helicon discharges, 
the inertia of electrons is found to strongly affect the high order radial modes. [5T] Actually, 
it has been verified experimentally that the inertia of electrons brings in measurable effect 
on the dispersion relation of waves whose frequencies are at least as low as 0.1|o; C e|.|S3 
For m e 7 ^ 0, Eq. (1.28) has two roots according to the quadratic formula: 


A 


k z + yjk 2 z - 4 6 rt kl 

25rt 


kz - \Az - 4 5 rt kl 
2 5 rt 




(1.39) 
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They can be simplified for 5 r tk -C k 2 : 


fa 


k 2 

"'W 

kz 


u ia 0 n e oq e 

'kz Bo ’ 


@2 


kz 

fat 


= -W c 

UJ 


(1.40) 


with Pi standing for the helicon mode (same as P, see Eq. ( 1 .29| ) ) and @2 evidently for 
an electron cyclotron mode (or TG mode, shortened for the initials of Trivelpiece and 
Gould who first discussed it with cylindrical boundaries[39]). The solution of Eq. (1.27) 
is eventually the combination of the following two Helmholtz equations: 


T"Bi a + /5{Bi a — 0, V“Bit) + ^2Bib — 0, 


(1.41) 


where Bi = Bi a -f Bib- A complete boundary condition for Eq. (1.41) was given by Chen 
and Arnush[33] for a conducting cylinder 


PiTiJ m (Tia) [(p2 + k z )Jm-i(T2a) + (p 2 ~ k z )J m+ i(T 2 a)] = 
faT2J m (T 2 a) [(Pi + k z )Jm-i{T 1 a) + (Pi - k z )J m+ i(Tia)] 


with T 2 = P\ — k 2 and T 2 = fa — 
and P 2 for a given mode structure, 


k 2 . From Eq. (1.42), we can solve the two roots of Pi 
and then obtain the corresponding dispersion relation. 


Helicon waves in a cylindrical non-uniform plasma 

In this section, we consider a plasma cylinder with radial density gradient, which is 
more realistic for helicon discharges than a uniform density profile that has been used 
so far. [19] The effect of the radial density gradient on helicon discharges was first recog¬ 
nised experimentally by Lehane and Thonemann[4C| and theoretically by Blevin and 
Christiansen [36] E], and examined in more detail by Chen and his colleagues. f42j [43] 
Although Blevin and Christiansen achieved an analytical dispersion relation for a specific 
radially non-uniform density profile, most previous studies rely on numerical computations 
for arbitrary radial non-uniformity in plasma density. In 2000, Breizman and Arefiev pre¬ 
sented a theoretical analysis to explore the effect of the radial density gradient on helicon 
waves, and discovered a new surface-type helicon mode. Though their analysis is based on 
a step-like radial profile of plasma density, the analytical dispersion relation they obtained 
is actually valid for general radially non-uniform density profiles, as long as the plasma 
column is long, thin (k z 1/a) and sufficiently dense (|w ce | <C u pe ). The surface-type 
helicon mode is called so because it is caused by a localised “surface” current, formed 
by the radial density jump. This radially localised helicon (RLH) mode may shed light 
on the high efficiency of helicon wave production of plasmas. Since the RLH mode will 
be employed in Chapter [4] to form a gap eigenmode in a plasma cylinder with periodic 
magnetic field, a short overview is given below. [35] 
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For \J LO C i | uJce | <w< | w ce | <C u; pe , we rewrite Eqs. (1.9)-(1.10) as: 


, ,2 , ,2 , ,2 

S = Hr> D = P = —f, 




ww r 


UT 


(1.43) 


and combine Eq. (1.1) with Eq. (1.11) to get the linear wave equation 

a; 2 

V x V x Ei = e • Ei. 


(1.44) 


The three components of Eq. (1.44) in cylindrical coordinates are: 


(™? , ,2^ c , im ( dE \0 , Ei e \ | , f _ cu 2 . nE , x ^ 

( 7T + k zJ E ir + — + lkz ~Q^~ ~ ~^2 ( SElr ~ iDE ie) , (1-45) 


• 9 f El A 9 ( dE ie Eio N 


l-) - y l -w-b —) + A: 2 Pi0 - k z —Ei z = ^ (S'Eie + iDEi r ) , (1.46) 

or \ r J or \ or r J r c A 


ik z d m 1 3 / <9Ei z \ m“ „ u; 


— — (rEi r ) - k z — E\o -— r- 

r or r r or \ or 


H- ^~E\ Z — -^-PE\ Z . 

r z c z 


(1.47) 


For a small longitudinal wave number (k z -C min(u; pe /c; 1/a)) and a sufficiently dense 
plasma (u -C |w ce | -C uj pe ), both Eq. ( 1.45|) and Eq. (1.46) read 


im 

r 

r 


or r 


(1.48) 


As will be shown, the radial nonuniformity of plasma density has a surprisingly strong effect 
on the structure of helicon modes with m / 0, and this effect is most pronounced in the 
limit of small k z considered here. For m ^ 0, expanding the quasineutrality (V ■ Di = 0) 
gives 

1 o i m 

- — [r (SE lr - iDE w )] + — (SE ie + iDE lr ) + ik z PE lz = 0. (1.49) 

r or r 


Equations ( 1.47 )-( [L49| ) are combined to form a closed set which can be reduced to: 


1 d I dE\ m 2 uo 2 

r dr V dr J r 2 c 2 lz ’ 


(1.50) 


1 d_ 

r dr 


Sr ^ (E lz - E) 


rn f rn S - 9E> \(E lz -E)-k 2 z PE lz = 0. (1.51) 

r V r dr J 


Here, a new unknown function has been introduced, E = E\ z — (k z r/m)Eig, which mea¬ 
sures the radial component of the perturbed magnetic field. Equations (1.50)- JL5T ) indi¬ 
cate two modes: a helicon mode with wave equation 


1 d ( dE\ m 2 

( r—— )- ttE = 


m u 


EdD/dr 


r dr \ dr J 


k 2 r c 2 1 — (mdD/dr)/(k^rP) 1 


(1.52) 
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and a TG mode with wave equation 


-|-(5r 
r or 


/IE 


lz 


dr 


m / m 


kiPE lz = 0. 


(1.53) 


For the helicon mode, Eq. (1.52) can be solved for u based on a step-like radial profile of 
plasma density 

n_, r < ro 
n+, r > r 0 


n, = 


(1.54) 


and the corresponding radial profile of electric field 


E — En 


(r/ro)l m l, r < ro 
(r/ro)“l” l l, r > r 0 


(1.55) 


where ro is the radius of the density discontinuity and Eq is a constant electric field. Under 
the limit of cu pe /c S> 1 /a, an analytical dispersion relation can be obtained 


m uj ce k 2 z c 2 


(1.56) 


The mode involves a perturbed “surface” current which is localised near the peak of the 
eigenfunction E{r) and distinguishes it from the modes studied previously. i'31, [41] This 
RLH mode can couple strongly to the antenna current during helicon discharges. 


For comparison, the dispersion relations of whistlers (Eq. (1.20)), helicon waves 
(Eq. (1.29)) and the RLH mode (Eq. (1.56)) are plotted together in Fig. 1.5 The plasma 



Figure 1.5: Dispersion relations of whistlers (uniform plasma in unbounded slab), helicon waves 
(uniform plasma in cylinder) and RLH mode (non-uniform plasma in cylinder) for a singly ionised 
argon plasma with peak density n s0 = 10 18 and field strength \B 0 \ = 0.01 T. 


parameters employed here are the same to those used for Fig. 0 For helicon waves, the 
root /? in Eq. (1.29) has been set to 3.83/a, taking the lowest root of Ji(Ta) = 0. For 
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the RLH mode, a numerical form factor T has been introduced to Eq. (1.56), which then 
has a new form w = Tuj ce c 2 k 2 /ujp e , for a Gaussian radial profile of plasma density. Fitting 
this new form with the computed dispersion relation results in T = —2.03.[55] Figure 1.5 


indicates that the dispersion curves of these three modes approach to each other when 
k z -» 0. 


1.1.2 Helicon discharge 


A helicon discharge usually refers to a plasma discharge driven by helicon waves which, 
as shown in Sec. [13 are RHCP electromagnetic waves propagating in a bounded magne¬ 
tised plasma with frequencies between the ion and electron cyclotron frequencies. mm 
Although helicon waves in a cylinder can be either RHCP or LHCP, the former is often 
dominant. mmmm A radio frequency (RF) antenna is usually employed to excite 
helicon waves and drive the helicon discharge. |49j Figure 1.6 shows the schematic of the 


first RF antenna (double-saddle) used to successfully drive a helicon discharge, [Tj together 
with those of other two widely used RF antennas: plane-polarised Nagoya Type iiiBiE] 
and twisted Nagoya Type III (half-turn helical). J3 : , [ 6 j The RF antenna is connected to a 



Figure 1.6: RF antennas for driving helicon discharges: (a) double-saddle antenna used to produce 
the first helicon plasma, [1] (b) plane-polarised Nagoya Type III antenna, mm, and (c) twisted 
Nagoya Type III antenna (half-turn helical). [H [5J '6j 

RF power supply and wrapped around an insulating chamber (usually a glass tube). The 
time-varying current on the antenna wire induces a time-varying magnetic field, which 
then induces an electrostatic field in the opposite direction of the current through Fara¬ 
day’s law. The induced electrostatic field is nulled by the build up of space charge. Legs of 
the antenna opposite to the machine axis have opposite induced space charge, giving rise 
to a transverse electric field. [48j This transverse field couples the transverse m = 1 mode 
structure, as shown in Fig. |1.3| (g\ = 1). The spatially oscillating electric field accelerates 
free electrons in the neutral gas and ionises it through collisional (high pressure) and col¬ 
lisionless (low pressure) heating mechanisms. [191 [50 j For an inductively coupled discharge 
which does not have a static magnetic field, the formed surface plasma shields the wave 
field from the antenna and allows the propagation of waves only within a short distance 
(known as skin depth) into the column. |T 2 l 49] If, however, an external magnetic field is 
applied axially through the plasma column, the dielectric scalar of the plasma becomes a 
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tensor (see Eqs. (1.8)-(1.10)), and waves of a certain frequency range will be allowed to 
propagate inside the plasma and not confined within the skin depth. mmm This fre¬ 
quency range of waves are typically helicon waves. Because they can propagate deep inside 
the plasma, which increases the efficiency of the wave energy deposition significantly, he¬ 
licon waves generate plasmas with densities typically an order higher than the inductively 
coupled discharge when operated at similar pressures and input RF powers. (32] 

An enhancement in collision frequency is usually required to fit experimental data to 
theoretical calculations for helicon discharges, jU 53, [53] and the high efficiency in helicon 
wave production of plasmas has not yet been completely understood. The pursuit of this 
understanding has stimulated a large amount of research and promoted the discovery of 
various applications using helicon plasma sources. mm These include, to name a few, 
processing semiconductor circuits, [32, 19] plasma propulsion for space travel, [55) 56, 5J, 58] 
gas laser media and plasma lenses for high energy particle beam, [32] magnetised plasma 
opening switches, [35] laser plasma sources, [59] and possibly electrodeless beam source and 
laser accelerator. [48] The helicon plasma can be also used as a research tool for measur¬ 
ing Hall coefficients in semiconductors. [60j studying Landau and Cherenkov damping, [60] 
measuring Doppler-shifted cyclotron damping. [61] studying magnetic fusion, [62] driving 
RF current, m, studying intense-beam plasma interactions, [35] understanding ionospheric 
phenomena, [601 [35] and studying Alfven wave propagations. |64] 

Helicon wave production of plasma was first achieved by Boswell in 1970,p] follow¬ 
ing ten years history of helicon propagation study first in solid state plasmas[651 [66 , 3D] 
and then in gaseous plasmas. m The basic theory of helicon waves was also studied ex¬ 
tensively in the 1960s. [M ESI EH EQl EH The early history of helicon studies prior to 
1997 is reviewed by Boswell and Chen. USUIS] Two questions remain unresolved: the rea¬ 
son for the high efficiency of helicon discharges and the dominance of the RHCP mode 
over the LHCP mode. [T9] This thesis first models the wave instability in helicon plas¬ 
mas with consideration of plasma flows, which have been observed recently in multiple 
helicon devices. EIE2] Second, the thesis studies the propagation of helicon waves inside 
a pinched plasma and supports the hypothesis that multiple radial modes could be ex¬ 
cited simultaneously, [T8l [731 ; 731 [75] which may contribute to the understanding of the 
“strangely” high efficiency. Due to the robust field line patterns of m = 1 mode, which 
does not vary with axial position as discussed in Sec. 1.1.1, helicon waves (particularly 
the RLH mode) will be utilised to form a gap eigenmode in the downstream region of a 
helicon discharge. 


1.2 Magnetic geometry 

The magnetic geometry of plasma confinement varies from uniform to focused, diverged, 
rippled, toroidal, dipolar, irregular et al., depending on the specific applications. Here, 
we consider four geometries: uniform as typically used for the helicon discharge, focused 
to get a pinched plasma for material processing, rippled to form a spectral gap and a 
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gap eigenmode through breaking the periodicity locally, and toroidal mainly for magnetic 
confinement fusion where the gap eigenmode found in the rippled field has applications. 


1.2.1 Uniform: WOMBAT 

For a uniform magnetic geometry, we consider the WOMBAT, [2D] which was designed to 
provide a plasma environment to study wave physics including chaos, turbulence, wave 
saturation and the related driving mechanisms. This device will be employed in Chapter [2] 
for investigating drift waves in a flowing plasma cylinder.[76] Figure [L7| shows a schematic 
and the employed magnetic field. It comprises a glass source tube and a large stainless 

0.04 p—-—-—-—-—i—-—-—-—-—i—-—-—-—-—i—-—-—-—-—i—-—-—-—-—^ 
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Figure 1.7: A schematic of the WOMBAT (Waves On Magnetised Beams And Turbulence) 
plasma experiment and the employed magnetic field. [7] 

steel diffusion chamber, which are attached to each other on axis. The source tube is 
50 cm long and 18 cm in diameter, while the diffusion chamber is 200 cm long and 90 cm 
in inner diameter. A large solenoid is employed inside the diffusion chamber to provide 
a steady magnetic held, which can be up to 0.02 T and highly uniform along the axis. 
To generate the plasma, a single loop antenna (20 cm in diameter, 1 cm wide and 0.3 cm 
thick) surrounding the source tube coupled the RF power to the argon gas. The 7.2 MHz 
electric current passing through the loop creates a time varying magnetic held, which in 
turn induces azimuthal current in the argon gas and leads to break down and formation of 
the plasma. [60] The base pressure of the diffusion chamber was maintained at 4 x 10~ 6 Torr 
by a turbomolecular pump and a rotary pump. In this work, the large solenoid was used 
to produce a high density blue core in the diffusion chamber. The dispersion relation of 


resistive drift waves, different from that of helicon waves shown in Eq. (1.20), Eq. (1.29) 
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and Eq. (1.56), is given in Eq. ( |2.12 ) coming from a two-fluid model. mm 

The PCEN (Plasma CENtrifuge[?9|). for which a two-fluid electrostatic flowing plasma 
model (used in Chapter [2]) was developed, m is a typical vacuum arc centrifuge designed 
for separating metal isotopes. It also employs a uniform magnetic field but with faster 
normalised rotation frequency, higher temperature and higher axial velocity than those of 
the WOMBAT. Moreover, the plasma source is not a helicon discharge but a DC discharge, 
namely the plasma is generated by a discharge between cathode plate and anode mesh. A 
schematic of the PCEN can be found in m 


1.2.2 Focused: MAGPIE 

To illustrate a focused magnetic geometry, the MAGPIE is employed. It is a linear plasma- 
material interaction machine which was recently built in the Plasma Research Laboratory 
at the Australian National University, and designed for studying basic plasma phenomena, 
testing materials in near-fusion plasma conditions, and developing potential diagnostics 
applicable for the edge regions of a fusion reactor. [8] Similar to other helicon devices, [46] 
MAGPIE mainly consists of a dielectric glass tube surrounded by an antenna, a vacuum 
pumping system, and a gas feeding system, together with a power supply system connected 


to the antenna, and various diagnostics. Figure 1.8 shows a schematic and the employed 
magnetic field, and introduces a cylindrical (r, 8, z ) coordinate system. The plasma is 
formed in the region under the antenna (—0.243 < z < —0.03 m) and the near field to the 
antenna. m Following convention, however, we define the whole glass tube (—1 < z < 0 m) 
as the source region and the compressed field region (0 < z < 0.7 nr) as the target region 
(or equivalently “diffusion region” in some references). In MAGPIE, the z < —0.243 m 
region is named “upstream” and z > —0.03 m “downstream”. A glass tube of length 1 m 
and radius 0.05 m is used to contain source plasmas in MAGPIE. A left hand half-turn 
helical antenna, 0.213 nr in length and 0.06 m in radius, is wrapped around the tube and 
connected to a tuning box which can be adjusted between 7 and 28 MHz, a directional 
coupler, a 5 kW RF amplifier, and a 150 W pre-amplifying unit. For the study presented 
in Chapter[3j a RF power of 2.1 kW, a frequency of 13.56 MHz, a pulse width of 1.5 ms and 
a duty circle of 1.5 % is used. The antenna current is measured by a Rogowski-coil-type 
current monitor. For these experiments an antenna current of magnitude I a = 38.8 A was 
measured. A grounded stainless steel cylindrical nresh surrounding the whole source region 
is employed to protect users. The source region is connected on-axis to the aluminium 
target chamber which is 0.7 m in length and 0.08 m in radius. Gases are fed through the 
downstream end of the target chamber, and drawn to the upstream end of the source tube 
by a 170 L/s turbo pump. Gas pressures are measured in the target chamber by a hot 
cathode Bayard-Alpert Ionization gauge (< 0.01 Pa), a Baratron pressure gauge (0.01- 
10 Pa) and a Convectron (0.1 Pa-101.33 kPa). In this experiment, argon gas is used with 
a filling pressure of Pb = 0.41 Pa. The two regions, source and target, are surrounded 
by a set of water cooled solenoids, with internal radius of 0.15 m. These source and 
target sets of solenoids are powered by two independent 1000 A, 20 V DC power supplies, 
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I (m) 


Figure 1.8: A schematic of the MAGPIE (MAGnetised Plasma Interaction Experiment), together 
with the employed magnetic field. 0 A circle denotes the position of the helicon antenna which is 
left hand half-turn helical. The dot-dashed line is the machine and coordinate system axis, defining 
r = 0 m. The coordinate system is right-handed with 9 = 0 chosen to be the zenith angle. The 
dotted coil pair is an optional coil to control the magnetic field gradient, which is not used in the 
present work. 


providing flexibility in the axial configuration of the static magnetic field, e. g. maximum 
of 0.09 T and 0.19 T in the source and target regions, respectively. The non-uniform 
field configuration is expected to provide a flexible degree of radial confinement, better 
plasma transport from the source tube to the target chamber, and possible increased 
plasma density according to previous studies. B m EH HI EDI EQ Wave propagations in 
the pinched plasma of the MAGPIE will be studied in Chapter [3| [54j 

1.2.3 Rippled: LAPD 

We use the LAPD to illustrate a rippled magnetic geometry. The LAPD is a large, linear 
plasma research device designed to study space plasma processes, |9] and more recently 
fusion research. It is ideal for studying the fundamental properties of plasmas on a large 
scale length, e. g. propagation of low frequency whistlers and SAW. wm As shown 
in Fig. |1.9[ the LAPD consists of four stainless-steel chamber sections, each of which is 
mounted on a wheeled platform to allow separate maintenance. The diameter and length 
of the machine are 1 m and 20.7 m, respectively. |10j The first bell-shaped section mounts 
a heater, a cathode and an anode. The heater is made of hollow alumina ceramic tubes 
which are threaded with tungsten filaments. This special design is expected to provide 
a uniform heating for the cathode with temperature slightly above the cathode emission 
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Figure 1.9: A schematic of the LAPD (LArge Plasma Device) and the employed magnetic field. [5] 
HO] 


temperature. When applied a negative voltage, the heated cathode emits electrons through 
both the thermionic electron emission and the DC discharge (towards the anode). The 
emitted electrons then go through the grid anode and ionise the neutral gas filled in 
the system, forming a plasma column along the system’s axis. The plasma column is 
18 m long and with densities up to 2 x 10 19 m~ 3 .(9j Radial expansions of the plasma are 
confined by an externally applied magnetic field, which is generated by a set of 68 magnet 


coils surrounding the chamber, as shown in Fig. 1.9 The held strength can be up to 
0.3 T, and varied with axial position by varying the current on each of the magnet coils, 
making it capable to configure a periodic array of magnetic mirrors. m Therefore, the 
LAPD is a promising device for identifying a spectral gap of plasma waves, and also an 
appropriate candidate for observing a gap eigenmode. [35] Discussions about identifying 
the gap eigenmodes of RLH waves and SAW will be given in Chapter [4] and Chapter [5j 
respectively. Moreover, there are 450 radial ports on the vacuum chamber thus providing 
excellent access for diagnostics. [9j 


1.2.4 Toroidal: tokamak 

Finally, a toroidal magnetic geometry is illustrated in term of a tokamak. Controlled ther¬ 
monuclear fusion can be an environmentally attractive and sustainable energy source, pro¬ 
viding a large scale base-load power. [82] However, realising the controlled thermonuclear 
fusion faces many challenges, which in science is the combined requirement of confining a 
sufficient quantity of plasma for a sufficiently long time at a sufficiently high temperature 
to make net fusion power. The Lawson criterion formulates this requirement specifically 
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m 


-3 


s keV, with h and T the peak ion density and temperature in the 


plasma and te the energy confinement time. [831184j We focus on the magnetic confinement 
fusion (MCF): the most achievable scheme. To avoid the loss of charged particles parallel 
to the magnetic field line, the magnetic field is usually shaped into a torus. A tokamak 
is an axisynunetric torus with a poloidal magnetic field that is produced by a toroidal 
electric current flowing inside the plasma, and a very strong longitudinal field parallel 
to the current. [55] The word “tokamak” is a transliteration from an acronym of Russian 
words: toroidalnaya kamera and magnitnaya katushka, meaning “toroidal chamber” and 
“magnetic coils”. [STj Figure 1.10 shows a schematic of the tokamak: two components of 
the magnetic field (poloidal and toroidal) forming a helical magnetic field nest confin¬ 
ing the plasma torus. The world’s largest tokamak is being under construction at the 



Figure 1.10: A schematic of the tokamak: poloidal (solid) and toroidal (dashed) magnetic field 
lines (a) form a helical magnetic field nest confining the plasma torus (b). 

Cadarache facility in the south of France: International Thermonuclear Experimental Re¬ 
actor (ITER), which aims to demonstrate the scientific and technological feasibility of 
producing commercial energy from the controlled thermonuclear fusion. It is designed to 
produce 500 MW of fusion power from 50 MW of input power, with the energy gain of 
Q = 10. [86] 

Given that fusion plasmas are energetically rich, complex exothermic physical systems, 
charged particles can still be lost even if the toroidal magnetic configuration is perfectly 
designed. Populations of charged particles can be accelerated well above thermal speeds 
through processes such as resonant wave absorption, charge exchange with injected high 
energy neutral beams, magnetic reconnection, and the fusion process itself. The impor¬ 
tance of these energetic populations is that they can drive electromagnetic waves which, 
in turn, can eject the same driving particles from confinement. m The most easily excited 
modes usually reside in their spectral gaps, which are formed by the various periodicities 
in the toroidally confined plasma. m These gap eigenmodes are introduced by symmetry¬ 
breaking due to toroidicity, plasma ellipticity and higher order shaping effects. Therefore, 
it is of practical interest to study the formation of gap eigenmodes and their interactions 
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with energetic particles. We shall investigate two types of gap eigenmodes: RLH waves in 
the whistler frequency range, which can be driven unstable by energetic electrons, [35] and 
SAW with frequency below the ion cyclotron frequency, which can be excited by energetic 
ions in tokamaks. 

1.3 Aims of the thesis 

The aim of thesis is to study the impact of magnetic geometry on spontaneous and driven 
wave modes in cylindrical helicon-source plasmas. Three questions will be addressed: 

(1) what are the plasma-driven wave modes in a rotating, flowing plasma with uniform 
magnetic field? 

(2) how do wave propagation characteristics of antenna-driven modes change with 
ramping magnetic field, and can this non-uniformity enhance the plasma density? 

(3) can gap eigenmodes be formed in a plasma cylinder with rippled magnetic field, 
by breaking the field periodicity locally? 
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Chapter 2 


Drift waves in a uniformly 
magnetised plasma with flows 


A two-fluid model developed originally to describe wave oscillations in the vacuum arc 
centrifuge, a cylindrical, rapidly rotating, low temperature and confined plasma column, 
is applied to interpret plasma oscillations in a RF generated linear magnetised plasma 
(WOMBAT) with similar density and field strength. Compared to typical centrifuge plas¬ 
mas, WOMBAT plasmas have slower normalised rotation frequency, lower temperature 
and lower axial velocity. Despite these differences, the two-fluid model provides a consis¬ 
tent description of the WOMBAT plasma configuration and yields qualitative agreement 
between measured and predicted wave oscillation frequencies with axial field strength. In 
addition, the radial profile of the density perturbation predicted by this model is con¬ 
sistent with the data. Parameter scans show that the dispersion curve is sensitive to 
the axial field strength and the electron temperature, and the dependence of oscillation 
frequency with electron temperature matches the experiment. These results consolidate 
earlier claims that the density and floating potential oscillations are a resistive drift mode, 
driven by the density gradient. To our knowledge, this is the first detailed physics model 
of flowing plasmas in the diffusion region away from the RF source. Possible extensions to 
the model, including temperature non-uniformity and magnetic field oscillations, are also 
discussed. 


2.1 Introduction 

Oscillations have been observed in the electric probe signals of many laboratory plasmas, 
including theta pinches,|88| [89J magnetic mirrors, [mum [ 92 ] Q machines, [93] vacuum arc 
centrifuge |£0| and helicon plasmas. [Ml [96UQlj 98 jj 99 ] Theta pinches and magnetic mirrors 
behave vastly different from the latter because of their high beta and significant magnetic 
curvature, which result in highly different plasma and field geometry. In contrast, Q 
machine, plasma centrifuge and helicon plasmas have low beta and negligible magnetic 
curvature, and have similar plasma properties. We exploit this similarity by deploying a 
two-fluid collisional model, initially developed to describe vacuum arc centrifuge plasmas, 
to describe the configuration and oscillations in a helicon plasma. Our goal is to better 
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Drift waves in a uniformly magnetised plasma with flows 


describe a class of electrostatic oscillations observed, but not yet conclusively identified, 
in the diffusion region of helicon plasmas, away from the RF source. 


The helicon plasma we study is provided by the WOMBAT device, whose schematic 
has been given in Fig. 1.7 The WOMBAT has been used to study helicon waves . [TOOi. 
E3 m] which have attracted great interest in the past decades due to their ability to 
produce high plasma densities. 11011 Such plasmas also provide a rich environment for the 
experimental study of drift waves and their impact on anomalous transport. [102 ] which is 
important to understand in fusion plasmas. Although there have been many publications 
on helicon wave physics, less attention has been devoted to the physics of low frequency 
oscillations in the diffusion region away from the source. |103 j Light et al. |10'41 [1051 
observed a low frequency electrostatic instability only above a critical magnetic field and 
identified it as a mixture of drift waves and a Kelvin-Helmholtz instability. Degeling et al. 
fOTHTM] observed relaxation oscillations in the kilohertz range that were associated with 
the various types of mode coupling in helicon discharges. More recently, a cylindrical linear 
model which treats global eigenmodes [ 107] has been employed to study drift waves in the 
magnetically confined plasma device VINETA. |1031 102 , 1108 ] however, the model does 
not involve the E x B rotation and neglects ion parallel motion which may influence the 
characteristics of low frequency oscillations. [ 109] Sun et al. [993 observed low frequency (13 
kHz) oscillations near the end of helicon source region, and identified them as resistive-drift 
Alfven waves. Sutherland et al. [110 ] observed low frequency ion cyclotron waves that were 
highly localised along the axial center of the WOMBAT plasma device. Their analysis of 
spectral measurements suggested the possible existence of a four wave interaction, where 
energy is down-converted to the ion cyclotron mode from the helicon pump. 


Low frequency oscillations in helicon plasmas can be broadly classified into two types: 
the Kelvin-Helmholtz instability and drift waves. The Kelvin-Helmholtz instability is 
driven by a velocity shear in a mass flow or a velocity difference across the interface 
between two fluids. m This can occur in helicon plasmas, because the streaming velocity 
of the ion fluid is much slower than that of electrons. The drift wave is a universal 
instability driven by a plasma pressure gradient perpendicular to the ambient field. [12j It 
can arise in fully ionised, magnetically confined and low-beta plasmas, m and has been 
observed in both linear and toroidal field geometries. jll3l fl 141 il 15l 1116.UT7 ] In WOMBAT 
plasmas, the velocity shear is small and the density and temperature gradients produce 
large pressure gradients, suggestive of large drive of drift waves. We thus restrict attention 
to this class of modes. 


In this chapter, a two-fluid model, which was developed originally for the vacuum arc 
centrifuge by Hole et al. [77, 118 ] to explain oscillations observed in the density and electric 
potential, is applied to study low frequency oscillations observed in the WOMBAT. We 
show that the equilibrium and perturbed density profiles and the space potential profile 
from the model and data are consistent. The model predicts unstable modes with a similar 
global mode radial structure, and a frequency comparable to observed signals. The chapter 


is organised as follows: Sec. 2.2 gives a brief description of the diagnostics on WOMBAT, 
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together with typically measured parameters, Sec. 2.3 introduces the two-fluid model and 


Sec. |2.4| discusses the wave physics revealed by this model and the data. Finally, Sec. |2.5 
discusses the possible extensions to this model and presents concluding remarks. 


2.2 Diagnostics and measured parameters 

To measure spatial and temporal profiles in the plasma, an uncompensated, translating 
Langmuir probe was inserted radially into the diffusion chamber, 50 cm from the source 
and diffusion chamber interface. The central wire of the probe was fed through an alumina 
support which was in turn shielded by a 6 mm diameter grounded steel tubing covering 
the whole extent of the probe length up to the probe tips. The removable probe tip was 
made of a 0.2 mm diameter and 8 mm long nickel wire. Radial translation of the probe 
was set using a computer controlled stepper motor arrangement that allowed the probe tip 
position to be selected with an accuracy of a few micrometers. To determine the plasma 
I(V) characteristics, the bias voltage on the Langmuir probe was swept between —50 V 
and +30 V using a Labview program. The plasma density (nj), electron temperature (T e ), 
plasma potential (V p ) and floating potential (Vf) were determined from the I(V) charac¬ 
teristics of the cylindrical probe. For the present high-density plasma discharge, the effect 
that RF fluctuations have on the floating potential, and hence the I(V) characteristics, 
was found to be negligible. The reproducibility for T e was within 0.5 eV, and the signal- 
to-noise ratio was improved by averaging many I(V) curves. The effect that fluctuations 
due to waves will have on time-averaged I(V) characteristics will be far more dominant 
than the RF fluctuations. These wave fluctuations would lead to error magnification in the 
second derivative of the I(V) curve and hence the electron energy distribution function. 
For this reason, the plasma parameters were obtained directly from the I(V) curve and 
not by the Druyvesteyn method. |119] 

Typical Langmuir probe measurements of WOMBAT argon plasmas are shown in 


Table 2.1 Other parameters, obtained for similar helicon plasmas. jl20 . 72] and measured 
by laser induced fluorescence, include: the ion temperature Tj. the bulk rotation frequency 
wo, and axial streaming velocity V z o- The parameter m, is the argon ion mass and the 
charge number Z is taken to match similar argon plasma conditions. [ 121 1 The remaining 


parameters are introduced in section 2.3 The antenna driving frequency is 7.2 MHz, 
which is much higher than the bulk rotation frequency and wave frequency. Finally, 
plasma parameters for the PCEN device are taken from published work. [791177] It should 
be noted that the operating gas in the PCEN work is magnesium while here it is argon. 


2.3 Plasma model 

2.3.1 Model assumptions 

The two-fluid model developed by Hole et al. HZ1I781 is based on the following plasma 
assumptions: 
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Table 2.1: Typical parameters of WOMBAT and PCEN (Plasma CENtrifuge). 


parameter 

WOMBAT 

PCEN [791 E2 

rii (on axis) 

1.4 x 10 19 nr -3 

5.2 x 10 19 m - 3 

T e 

1.5 eV 

2.9 eV 

Ti 

0.1 eV 

2.9 eV 

mi 

40 amu (Ar) 

24.31 amu (Mg) 

B z 

0.0185 T 

0.05 T 

Z 

1.0 

1.5 

V z0 

200 m s~ x 

10 4 m s _1 

U 0 

1.42 krad s _1 

184 krad s _1 

, , _ B z eZ 

~ m . 

44.5 krad s ~ 3 

295 krad s _1 

n i0 = ^ 

0.032 

0.59 

= (3 +z) kB ? e D2 

v T e ' miUJ^Rr 

0.54 

1.6 

X — eZriio vl 

°rs ~ Bz 7B 

0.0235 

0.03 

^(characteristic radius) 

6 cm 

1.43 cm 


(1) Ions of different charge can be treated as a single species with average charge Z. 

(2) The plasma is quasi-neutral, giving that n e = Zn t . 

(3) The steady-state plasma is azimuthally symmetric and has no axial structure. 

(4) The effects induced by plasma fluctuations on the externally applied field can be 
neglected. 

(5) Both finite Larmor radius (FLR) and viscosity effects are negligible. 

( 6 ) For the range of frequencies considered here, the electron inertia can be neglected. 

(7) The ion and electron temperatures, T t and T e , are uniform across the plasma 
column. 

( 8 ) The steady-state ion density distribution has a form of no = n*(0)e^ r /' R ) 2 , which 
is a Gaussian profile. Here, nj(0) is the on-axis ion density, and R is the characteristic 
radius at which the density is 1 /e of its on-axis value. 

(9) The steady-state velocities of ions and electrons can be written as Vi = (0,u^r, Uj 2 ) 
and v e = ( 0 , u e (r)r, u e2 (r)), respectively, where w* is the ion rigid rotor rotation frequency, 
Vi z is the ion uniform axial streaming velocity, oj e (r) is the electron rotation frequency, and 
v ez {r) is the electron streaming velocity. While treated in other work [78] . radial diffusion 
of both ions and electrons due to electron-ion collision is negligible. 

Here, length and time are normalised to R and 1 /u>d respectively. A normalised 
cylindrical coordinate system is then developed, with (?”, 0, z) = ( r/R , 6 , z/R) and 
t = uj c it, where f and z are the normalised radial and axial positions, respectively. 
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2.3.2 Two-fluid equations 

The model comprises the motion and continuity equations of ion and electron fluids, 
written respectively as 


chii 

Ih 


+ (ui • V)ui = + ArVIj) + u ; x z + 5 rs n c iC • (u e - Ui), 


ipZ(-S7li + Vx) - u e x z + A rs n cZ £ • - u e ) = 0, 

dk „ 

- — = V • ui + ui • VZ;, 
dk 


- — = V • U e + U e • Vli, 
or 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 


with terms defined as follows: 


U; = —^ = {r<pri, rfli,Uiz), u e = 


XV'7 


OJrjR 


— ( r Prei U ez ), 


At = w, V = 

J e rriiUJ, 


kBT e 


,X = 


e 0 


k = In 


2 « 2 ’ 

n,; Z7 




= diag(^_L,C±,l), 


nj( 0 ) 


5 rid 


^ei(O) 




eZriio rj L 

B z 'Ye 


Here, subscript z and e refer to ion and electron parameters respectively, is the nor¬ 
malised radial velocity divided by f, H is the normalised rotation frequency, u z indicates 
the normalised axial velocity, A t is the ratio between ion and electron temperatures, ip is 
a convenient constant which for At = 1 becomes the square of the normalised ion thermal 
velocity, \ is a normalised electric potential (f>, k is the logarithm of the ratio of the ion 
density rq to its on-axis value ni(0), and n c i is the ratio of the electron-ion collision fre¬ 
quency v e i to its on-axis value 1 ^( 0 ). Also, S rs is the normalised resistivity parallel to the 
magnetic field, where tjl is the electrical resistivity of a Lorentz gas and 7 e is the ratio of 
the conductivity of a charge state Z to that in a Lorentz gas. [122 j 


2.3.3 Steady-state solution 


In cylindrical geometry with a purely axial constant field B = (0,0 ,B Z ), the steady-state 
solution of this model is given by 


Xo(r) 


Xc + 


Hjo 
2 ipZ 


(1 + Hjo) + 


A t 
Z 


f 2 . 


(2.5) 


with 

H e o = t2jo(l + Hjo) + + Z ), (2-6) 

where Xc is an arbitrary reference potential. The axial current in this model is uncon¬ 
strained, and can arbitrarily be set to zero (ui Z o = UezO ), consistent with WOMBAT 
boundary conditions. 
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We have compared the steady-state solution to experimental data. Measurements 
of the fluctuation frequency from probe measurements suggest rigid rotation, which is 
consistent with the spectroscopic data obtained by Scime et al. [72] Figure [2T| a) shows the 
equilibrium density profile of the WOMBAT plasma for two different axial held strengths, 
B z = 0.0185 T and B z = 0.0034 T, respectively. [IT] Overlaid are best fits to the profile 



6 

5 

4 

% 3 
2 

1 




• fl z =0.0185 T 
■ B,=0.0034 T 


(c) 


0.05 


0.10 

r(m) 


0.15 


0.20 


Figure 2.1: Plasma configuration of WOMBAT: (a) Gaussian fit of the equilibrium density profile, 
for B z = 0.0185 T and B z = 0.0034 T respectively: fill (b) the radial variation of the normalised 
space potential (the measurable electric potential is about = 8.6 x 10 _5 T e Xo V) from the fit 
procedure (solid line), the model (dashed) and the data (dots),[7] which was taken under the 
conditions of B z = 0.0068 T, 0.6 mTorr and 1708 W in the experiment; (c) the radial profile of the 
electron temperature at two field strengths: B z = 0.0185 T and B z = 0.0034 T. m The dashed 
vertical line in (c) corresponds to the radial location of fluctuation measurements. 


using a Gaussian profile (mo = m(0)e~^ r ^ R ^ 2 ): these show reasonable agreement over 
the range of the data. The stronger held produces better conhnement with a smaller 
characteristic radius R. Figure [27T| (b) shows the radial variation of the space potential in 
WOMBAT for a slightly different held strength, B z = 0.0068 T. [7] Overlaid is a best ht 
to the data using the parabolic potential prohle in Eq. (2.5), with the arbitrary reference 
potential (x c ) and gradient G = |^(1 + flio) + ^ free parameters. Although the ht, for 
which G = 0.5, is reasonable, the scatter in the data is large. For the same Xc, we have 
over plotted the model potential prohle with G = 0.04. Within the bulk of the plasma, 
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out to the characteristic radius r = 0.08 m, model and observed potential profiles broadly 
agree. Figure |2.1[ c) shows the radial profile of T e in the WOMBAT plasma. m The 
constant electron temperature of the model is an approximation to the radially varying 
experimental profile. For the model, we have chosen the T e value at the position at which 
the frequency was determined, r = 7 cm. 


2.3.4 Normal mode analysis 


To compute the normal modes of the system, Hole et al. Em apply a linear perturbation 
treatment with plasma parameters C taking the form 


C(r, f, 6 , z) = Co (f) + eCi (f)e^ me+ksS ~ WT \ 


(2.7) 


Here, e is the perturbation parameter, m the azimuthal mode number, k~ the axial wave 
number and u the angular frequency. We consider a complex value of ui so as to study 
both the mode frequency and growth rate, i. e. u = cu r + iui 1 . To first order in e the 


system of Eq. (2.1)-(2.4) reduces to 


( #'i (y)-X[(y)\ \ 

white) 

yp'feAv) - im^X[(y) 

y6 rs Z ± e-yX[(y) 

V 0 


= A 


( ki(y) \ 
Xi(y) 
<Pril(y) 
<Prel{y) 

V Uez(y) J 


( 2 . 8 ) 


where A is the matrix 

( 

A = 




m^/C 


0 

ivj _ iC 2 

iC 

0 

\ 



2tuy 


2 2zv 

2vj 

1 

2 

VJ 

_ ^ (m 2 
zu\ y 

+ fef) 

0 

-i + y-lf 

m 

2vj 

0 


K 

w 

~ - 

- 2m+ 

o 

0 

-1 + 2/ 

ikz 

2 




0 


2 y 

0 

1 

2 

0 




0 


—ikz'i 

0 

0 

0 

/ 


+ 6 r 


/ 

0 

0 

Ue~ y 

2 

Ue~ y 

2 

0 

\ 


0 

0 

0 

0 

0 



0 

0 

im£ ±e~ y 

2 

im£ ±e~ y 

2 

0 



Cue-n-f^ + ^f + +) 

0 

iCZ ± e~ y 

2vj 

+ e-» 
2vj 

0 


V 

e~ y k~'& 

VJ 

0 

0 

0 

5J 

1 

QJ 

1 

) 


(2.9) 


Here, w = oj — 77+2*0 — kzUzo is the frequency in the frame of the ion fluid, (7=1 + 212*0, 
+ = (Ay + Z)'ip, y = f 2 and we have introduced a new dependent variable X\ (y). where 


Xi (y) = 


z 


Z + A t 


[hi{y) - xi(y)]- 


( 2 . 10 ) 

























30 


Drift waves in a uniformly magnetised plasma with flows 


For large axial wavelength modes of the resistive plasma column, for which k 2 < 5 rs , this 
model can be reduced to a second order differential equation 


( m 2 -C 2 


L(N c )[ 9l (y)\ = 0, 


( 2 . 11 ) 


where 


L{Nc) - y h 2 +{l ~ y) h + 


N r . 



N c = 


(m 2 - C 2 )[m + §/(y)] , mC 


+ 


m 


m — mD 2 0 + i^f{y) 

and f(y ) = F 2 e y with the normalised axial wave number F = kz/\/S^. For odd m modes, 
the boundary conditions for Eq. (2.11) are gi(0) = 0 and gi(Y) = 0 with the infinite 
radius Y representing the edge of plasma column. For even m, these conditions become 
g\ (0) = 0 and gi(Y) = 0. We seek unstable solutions to Eq. (2.11), for which m l > 0. The 
solutions m = ±C are stable and hence discarded. 


2.4 Wave physics 

The experimental results indicate that a m = 1 drift wave mode dominates the density 
fluctuations, and that the mode appears to propagate purely azimuthally in direction 
of the electron diamagnetic drift. m There is no clear radial component of the mode 
propagation. The observed phase velocity of the mode is approximately 200 m/s in the 
region of maximum density perturbation. 


2.4.1 Dispersion curve 


As in Hole et al, m we have solved Eq. ( |2.11[ ) by a shooting method but using our new 
code written in Mathematica. For m = 1 the boundary conditions are g\ (0) = gi(Y) = 0. 
As the differential equation is homogeneous, the gradient at the edge g\ (Y) is arbitrary: 
we have chosen g[(Y) = 1. To solve for given F, we choose a trial m and march the 
solution from the edge to the core. The complex frequency m is then adjusted until the 
on-axis boundary condition is satisfied. We commence the procedure at F = 0, for which 
an analytical solution for m is available. This is found from the solution to 


m(m 2 — C 2 ) mC 
N r = — -H-= 2n + \m\, 


m — mP| 0 


m 


( 2 . 12 ) 


where n is the number of radial nodes in the plasma. Figure |2.2| shows the dispersion 


curve for the plasma conditions of the PCEN in Table 2.1 We have compared the disper¬ 


sion curve to Hole et al, m and found it to be identical, thus validating our numerics. 
For PCEN plasmas, the peak of normalised growth rate m l = 0.39 lies at F = 0.3, with 
normalised frequency m r = 0.007: the wave is thus near stationary in the frame of the 
ion fluid. The mode crossings located at F = 0.55 are associated with the centrifugal 
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Figure 2.2: Dispersion curves of PCEN (dashed line) and WOMBAT (solid line), generated 
by the model based on the conditions shown in Table 
measurable growth rate is about 2.95 x 10 5 cC rad s -1 


2.1 


(a) normalised growth rate vj 1 (the 
for PCEN and 4.45 x 10 4, cC rad s _1 for 
WOMBAT); (b) normalised frequency w r (the measurable frequency is about 2.95 x 10 5 (t*7 r + 
0.59m + 3.4 x 10~ 2 k z ) rad s _1 for PCEN and 4.45 x 10 4 (tA7 r + 0.032?n + 4.5 x 10 _3 fc z ) rad s _1 for 
WOMBAT) vs normalised axial wavenumber F = k~/\/5 rs (the measurable axial wavenumber is 
about k z = 12.IF m _1 for PCEN and k z = 2.6F m _1 for WOMBAT). 


instability. Also shown is the dispersion curve for WOMBAT plasma conditions of Ta¬ 


ble 2T For WOMBAT plasmas, the peak growth rate of w l = 0.14 occurs at F = 0.37, 
for which w r = 0.44. There is no mode crossing in the dispersion curve for WOMBAT 
plasmas. Also, w r > 0, and so in the laboratory frame of = vo r + mSl t o + k z u z o, approx¬ 
imately u r ~ vj t + fljO) such that the frequency of unstable waves is always larger than 
the sum of the plasma rotation frequency and axial velocity, and the wave propagates in 
the + 9 direction (electron diamagnetic drift direction). We attribute the difference be¬ 
tween dispersion curves to the very low ion temperature Tj and slow normalised rotation 
frequency fljo in WOMBAT plasmas, compared with PCEN plasmas. Indeed, this can be 
seen by varying the parameters from PCEN to WOMBAT conditions. If *P is held fixed 
and lowered, the growth rate drops linearly with flto, and the dispersion curve shifts 
to lower kz- Conversely, if D ? ;o is held fixed and T lowered, the dispersion curve expands 
to higher k z - and the growth rate drops. Thus, the decrease in growth rate is due to both 
the decrease in rotation and temperature, while the expansion of the dispersion curve to 
higher k z is due to the lower temperature. The corresponding increase in w r at low k z 
occurs because as rotation drops, the Doppler shifted frequency, w r = u r — — kgUg o, 

increases. 


2.4.2 Dispersion curve sensitivity with plasma parameters 

We have examined the sensitivity of the WOMBAT plasma dispersion curve with Djo, B z 


and T e . For Djo < 0.05 the dispersion curve changes by less than 1%. Figure T3 shows 
the change in dispersion curve for different B z and T e . Figure 


2.3 


a), (b) shows that w l 

and vo r decrease with increasing B z : the physical growth rate w* = w l u; c i also obeys this 
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Figure 2.3: Results of dispersion curve sensitivity with B z and T e : (a), (b) scanning results of B z - 
(c), (d) scanning results of T e . Here, the measurable growth rate is about 4.45 x 10 i w l rad s _1 and 
the measurable frequency is about 4.45 x 10 4 (w r + 0.032m + 4.5 x 10 ~ 3 k z ) rad s _1 , with the 
measurable axial wavenumber about k z = 2.6 F m -1 . 


trend while u r = ’uj r u c i exhibits the same trend in the range F > 0.5. Also, the peak 
growth rate shifts to lower k z , and so for constant R the axial wavelength of the most 
unstable mode is increased. In contrast, Fig. 2.3 c), (d) shows that vj 1 and -w r increase 
with increasing T e , and the peak growth rate shifts to larger k z . The dependence of growth 
rate with pressure gradient, and frequency with inverse field strength is consistent with a 
resistive drift wave. [12] 


2.4.3 Wave oscillation frequency with B z and T e 

In this section we compute the variation in wave frequency at the maximum growth rate 
with field strength, for two choices of electron temperature: a constant T e = 2.13 eV, 
independent of B z \ and a B z dependent T e . As T e data is only available at two field 
strengths (Fig. |2.1| (c)): B z = 0.0185 T, T e = 3.2 eV; and B z = 0.0034 T, T e = 1.5 eV, 
we have used linear interpolation to compute T e at other field strengths. Fluctuation 
data is available at five field strengths: B z = 0.0129 T, B z = 0.0146 T, B z = 0.0162 T, 
B z = 0.0179 T and B z = 0.0195 T. We have solved the dispersion curve for each case, and 
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selected w r corresponding to the maximum w l to calculate the frequency in the laboratory 
frame. The maximum error in the electron temperature measurements is ±0.5 eV which 
was determined from the reproducibility of the experimental measurements. This is also 
the minimum temperature that we could measure. 


Figure 2.4 shows the measured and predicted oscillation frequencies. For both curves, 



0 -■-■-■- 1 -■-■-■- 1 -■-■-■- 1 -■-■-■- 

0.012 0.014 0.016 0.018 0.020 

B z ( T) 


Figure 2.4: Wave oscillation frequency with B z and T e 


the predicted and measured frequencies differ by a factor of 3.5. There are various possible 
reasons for this gap. The most likely is that the model assumes flat T e and Tj profiles, 
whereas WOMBAT plasmas exhibit a non-uniform T e profile. It may also be the case that 
the Tj profile is non-uniform. Second, there may be slip between the rotation of WOMBAT 
plasmas and the oscillation frequency. Third, we have neglected the fluctuations in the 
magnetic field, which may affect the dispersion curve. Finally, the calculation of T e for 
each B z is a linear interpretation between the two known data, which may bring errors 


and imprecisions into the prediction. Figure 2.4 also reveals that two trends predicted by 
using a constant T e = 2.13 eV and a B z dependent T e are divergent. Predictions for the 
B z dependent T e profile exhibit the same trend as the data, suggesting the model may 
have correctly captured the T e dependence with B z . Finally, to ensure that this model 
is in principle capable to reproduce the data, we adjusted T e to fit the data points. We 
found that a value of T e = 0.5 eV at low field dropping to T e = 0.3 eV at high field was 
able to reproduce the observed frequency. 


2.4.4 Perturbed density profile 

Figure [±5] shows the measured and predicted radial profiles of the perturbed density n t \, 
together with the equilibrium density gradient \n' i0 (r)\. The data used here were taken 
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r(m) 


Figure 2.5: Comparison between the measured (dots) and predicted (solid line) radial profiles of 
the perturbed density Tin, together with the equilibrium density gradient |n' 0 (r)| (dashed line). 


under the conditions of B z = 0.01 T , R = 7.6 cm and T e = 2.46 eV. The linear perturbation 
treatment gives no solution for the absolute magnitude of nn , and so we have fitted for the 


amplitude. Inspection of Fig. |2.5| suggests m i has a single peak, and therefore consistent 
with n = 0 of the perturbed mode. Also, the peak of the eigenfunction occurs in the region 
where the equilibrium density gradient is large, suggesting that the instability has resistive 
drift-type characteristics. m The exact radial position of the peak in nn however does 
not agree. This could be due to uncertainty in the core position of the plasma, or the 
assumption of constant T e in the model. 


Finally, to complement our visualisation of the mode structure, we have computed 
the vector field of the linear perturbed mass flow by mi(nn\ii 0 + n,;oU;i). The perturbed 
velocity components u,;i = (ripfn, rflu, u^i), u e i = (ripfei, 7f2 e i, itezi) and perturbed 
density n t \ can be computed from the solution of g\(y) and following equations: 


ki (y) 


-5i (y) 

1 -1- i 

m^ 0 +2m^-ro 

1 ' 

f{y)-2im 


xi (y) = - (i + ■y)s'i(y), 
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^ii(y) = — \m +i[<^ fel (y) - CVfji(y)]} , (2.17) 

l y J 


Uizi(y) = '^~ bF ^ f In (y). 


VO 


(2.18) 


Figure 2.6 shows the flow vector field at time t = 0. As time advances, the mass flow 
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Figure 2.6: Vector field plot of the real part of the perturbed mass flow for mode (to = 1, n = 0) 
in a cross-section of the WOMBAT plasma. 

vector field rotates in the clockwise direction, which is the direction of electron diamagnetic 
drift with B z into the page. Here, the coordinate x and y label the cross-section of the 
plasma, with x = r cos 6 and y = rsin0. The 9 vector points clockwise with z into 
the page. The mass flow is zero on axis (x = 0, y = 0), consistent with the boundary 
conditions for density. The symmetric mass flow represents a flow around a stagnant core, 
with the largest flow in regions of peak density gradient. While we have only studied 
linear oscillatory response in this work, a similar nonlinear flow pattern may contribute 
to rotation damping. 

2.5 Conclusions 

In this chapter, we employed a two-fluid model, which was developed originally for de¬ 
scribing the wave oscillations observed in PCEN, [77| 78, 118 ] to study the low frequency 
oscillations in WOMBAT. To ensure that this model is consistent with WOMBAT, the 
measured and predicted plasma configurations were compared, including the equilibrium 
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density profile and the space potential profile. These show that although the density and 
space potential profiles agree, the temperature profile is not flat, as assumed by the model. 
Next, dispersion curves were generated for WOMBAT plasmas. Compared to the more 
rapidly rotating centrifuge plasma, the drift wave instability, unstable at larger wave¬ 
lengths, has a normalised frequency much larger than the normalised frequency for the 
centrifuge. The difference between dispersion curves is principally due to the low rotation 
speed of WOMBAT plasmas compared to those of the centrifuge. 

Our study of the wave oscillation frequency with B z and T e reveals that the mea¬ 
sured fluctuation frequency is a factor of 3.5 lower than the predicted frequency, and the 
predicted trend of oscillation frequency with B z for inferred T e matches the data. The dis¬ 
crepancy between the measured and predicted fluctuation frequencies may be attributable 
to the limitations of assumed temperature uniformity across plasma column. Weaker 
model limitations may be our neglect of effects induced by plasma fluctuations on the ex¬ 
ternally applied field. These effects are included in a TwO-fluid Electromagnetic Flowing 
pLasma (TOEFL) model developed in Appendix [Aj for which, however, the computation 
has not yet been finished. Data limitations include the simultaneous measurements of 
plasma rotation profile and probe frequency, and the measured dependency of T e with B z . 

Finally, we find that the measured and predicted perturbed density profiles have a 
single peak in the radial direction, indicating the perturbed mode has n = 0, and the 
peak position is in the region of maximum equilibrium density gradient. This qualitative 
agreement consolidates earlier claims that the mode is a drift mode, driven by the density 
gradient of the plasma. To our knowledge, this is the first detailed physics model of flowing 
plasmas in the diffusion region away from the RF source. 



Chapter 3 


Wave propagations in a pinched 
plasma 


A RF field solver j 123] based on Maxwell’s equations and a cold plasma dielectric tensor 
is employed to describe wave phenomena observed in a cylindrical non-uniform helicon 
discharge. The experiment is carried out on a recently built linear plasma-material inter¬ 


action machine: the MAGPIE (see Fig. 1.8),[SI in which both plasma density and static 
magnetic field are functions of axial position. The field strength increases by a factor of 15 
from source to target plate, and the plasma density and electron temperature are radially 
non-uniform. With an enhancement factor of 9.5 to the electron-ion Coulomb collision 
frequency, a 12% reduction in the antenna radius, and the same other conditions as em¬ 
ployed in the experiment, the solver produces axial and radial profiles of wave amplitude 
and phase that are consistent with measurements. A numerical study on the effects of 
axial gradient in plasma density and static magnetic field on wave propagations is per¬ 
formed, revealing that the helicon wave has weaker attenuation away from the antenna in 
a focused field compared to a uniform field. This may be consistent with observations of 
increased ionization efficiency and plasma production in a non-uniform field. We find that 
the relationship between plasma density, static magnetic field strength and axial wave¬ 
length agrees well with a simple theory developed previously. A numerical scan of the 
enhancement factor to the electron-ion Coulomb collision frequency from 1 to 15 shows 
that the wave amplitude is lowered and the power deposited into the core plasma decreases 
as the enhancement factor increases, possibly due to the stronger edge heating for higher 
collision frequencies. 


3.1 Introduction 

To date, most helicon studies have treated devices with uniform static magnetic fields, 
however, many applications require operation with axial magnetic field variations. m a 
few researchers have investigated helicon plasma sources with non-uniform magnetic fields, 
and have found that the plasma density increased when a cusp or non-uniform magnetic 
field was placed in the vicinity of the helicon antenna. [18, 07J [19, 80} However, detailed 
examination of the reasons for this enhanced plasma density has not yet been conducted, 
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although fast electrons and improved confinement are mentioned as possible contributors. 
Guo et al.[81j furthered this study by looking at the effects of non-uniform magnetic 
held on source operations, and found that a strong axial gradient in density associated 
with non-uniform held configuration can contribute to the absorption of wave helds and 
a high ionization efficiency. Takechi et al. j!25 j also suggested that there may be a close 
relationship between plasma density profile and RF wave propagation and absorption 
regions, hnding the density uniformity in the radial direction improved markedly with the 
cusp held. Therefore, studying the effects of various static magnetic held configurations on 
helicon wave propagation is of signihcant importance to producing desired plasma prohles 
and understanding the role of magnetic held in helicon plasma generations. 

This chapter is dedicated to modelling the wave held observed in MAGPIE, and inves¬ 
tigating helicon wave propagation in the non-uniform magnetised plasma of this machine, 
in which both the static magnetic held and its associated plasma density are functions 
of axial position. The plasma density and electron temperature are also dependent on 
radius. Independent measurements of electron temperature T e in MAGPIE at lower held 
and power conditions show that T e does not change substantially along We have thus 
assumed T e is independent of z in the present study. The static magnetic held is almost 


independent of r (Eq. (3.5)). A RF held solver,[124] based on Maxwell’s equations and a 
cold plasma dielectric tensor, is employed in this study. The motivations of our work are 
to explain the wave held measurements in MAGPIE, and to study the effects of magnetic 
held conhguration on helicon wave propagation. The rest of the chapter is organised as 
follows: Sec. |3.2| describes the diagnostic tools, together with the measured static mag¬ 


netic held, plasma density and temperature prohles; Sec. 3.3 provides an overview of the 
employed theoretical model and the numerical code, together with comparisons between 


computed and measured wave helds; Sec. 3.4 is dedicated to a numerical study of the 


effects of plasma density and static magnetic held prohles on the wave propagation char¬ 


acteristics; and Sec. 3.5 studies the physics meaning of the enhancement factor to the 
electron-ion Coulomb collision frequency, and the effects of the direction of static mag¬ 


netic held on wave propagations. Finally, Sec. 3.6 presents concluding remarks and future 
work for continuing research. 


3.2 Experiment 


3.2.1 Plasma profile diagnostics 

A passively compensated Langmuir probe was employed in our experiment to measure 
the plasma density and electron temperature, calculated from the I(V) curve obtained 
by an Impedans Data Acquisition system. [126] The probe comprises a platinum wire of 
diameter 0.1 mm, and a surrounding alumina insulator. The length of the insulator is 
6 mm shorter than that of the platinum wire so that the exposed platinum wire forms 
the probe tip. Electron currents were drawn to clean the probe during regular intervals of 


argon discharges. The probe is located at z = 0.17 m as shown in Fig. 1.8 
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Figure 3.1: Typical measured profiles: (a) axial profile of static magnetic field on axis, (b) radial 
profiles of plasma density (dots) and electron temperature (squares) at z = 0.17 m, together with 
their fitted lines, solid and dashed, respectively. The solid bar in (a) denotes the antenna location. 


Typical measured axial profile of field strength and radial profiles of plasma density 


and electron temperature in MAGPIE are shown in Fig. 3.1 As shown in Fig. |3.1[ a), 
the increase in field strength, Bq(z), from antenna end (z = —0.243 m) to field peak 
(z = 0.51 nr) is a factor of 15. The axial profile of plasma density, n e (z), is assumed to 
be proportional to Bq(z), consistent with generally accepted knowledge that the density 
follows the magnetic held linearly. [49, 127j Figure [3T| (b) shows the radial profiles of plasma 
density, n e (r), and electron temperature, T e (r), measured at z = 0.17 m and fitted with 
straight lines. During the density fitting procedure, in order to avoid negative fitted values, 
the density was set to zero in the region of 0.066 < r < 0.08 m. We assume the total 
density profile is separable, such that n e (r,z ) = n e (r) x n e (z). The fitted lines in n e (r) 
and T e (r), and the measured Bq(z) data will be used in Sec. 3.3 to constrain wave held 
simulations. 


3.2.2 Wave field diagnostics 

Helicon wave helds were measured by a two-axis “B dot” or Mirnov probe. Details about 
the probe can be found in Blackwell et ah. [8] To measure the axial prohles of B r and 
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B z , the probe was inserted on axis from the end of the target chamber. The probe is 
long enough to measure B r and B z in the range —0.25 < z < 0.7 m. Two perpendicular 
magnetic field components ( B r and B z in this case) can be sampled simultaneously. To 
measure the radial profiles of the three magnetic wave components, B r , Bg and B Zl the 
probe was inserted radially at z = 0.17 m, and rotated about its axis to measure Bg and B z . 
The B-dot probe couples inductively to the magnetic components of the helicon wave and 
electrostatically to the RF time varying plasma potential. To limit our measurements to 
the inductively coupled response, a current balun was employed to screen the electrostatic 
response. Further information about the procedure to eliminate the electrostatic response 
of the probe can be found in Franck et al.. m Both axial and radial profiles of wave 
phase were measured through a phase-comparison method, similar to Light et al.. [75J To 
measure the variation in wave phase with axial position, the signal from an on-axis axially 
inserted probe was compared to the phase of the antenna current. A similar procedure 
was conducted to measure the variation in wave phase with radial position at z = 0.17 nr. 
It should be noted that all probe diagnostics are intrusive, and can affect the plasma 
parameters, and hence the wave fields. 


3.3 Simulation 

A RF field solver (or ElectroMagnetic Solver, EMS |124) f based on Maxwell’s equations 
and a cold-plasma dielectric tensor is employed in this study to interpret the RF waves 
measured in MAGPIE. This solver has been used successfully in explaining wave phenom¬ 
ena in two other machines: a helicon discharge machine at The University of Texas at 
Austin [3Hj and the LAPD at the University of California at Los Angles. m Since these 
two machines have similar plasma density, temperature and particle species to those on 
the MAGPIE, the EMS should be applicable to describing the wave observations in MAG¬ 
PIE. Details of the solver can be found in Chen et al..[ T24 j while a brief overview is given 
below. 


3.3.1 EMS 


The Maxwell’s equations that this solver employs to determine the RF wave field in a 


helicon discharge are eventually Eq. (1.1) and Eq. (1.11) but including an external antenna 
current j a : 

V xE = iwB, (3.1) 


— V xB = -LuD + j E 
ho 


(3.2) 


Here, we have dropped the subscript “1” for simplicity. Because the wave field is driven 
by the antenna, the wave frequency is the same to the antenna driving frequency u>. The 


quantities D and E are linked via the cold-plasma dielectric tensor e (Eq. (1.8)) that 
represents vacuum, glass and plasma. In the vacuum and glass regions, the dielectric 
tensor is e = e*<5, where 5 is the Kronecker symbol and e* is a scalar. The term e* equals 
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to 1 and e g for vacuum and glass regions, respectively, where e g is the dielectric constant 
of glass. In the plasma region, the relation between D and E has the form [21] 

D = e 0 [SE - iD{ E x b) + (P - S)( E • b)b], (3.3) 


where b = Bo/Po is the unit vector along the static magnetic field. To take into account 
the collisions of particles in the discharge, we introduce a phenomenological collision term 


v s for each species into e, and rewrite Eq. (1.9) and Eq. (1.10) to: 


5 = 1-E 


ul s {u+iv a ) 

u[(uj+iv a ) 2 ~u} 2 s 


D = E 




u)\(u+iv a y 2 -CJ% S 


(3.4) 


p= 1-E 




c j)(uj-\-iv s )' 


For plasma parameters typical of MAGPIE, the electron-neutral collision frequency is 
found to be an order of magnitude smaller than the electron-ion collision frequency, and two 
orders smaller than the electron-ion collision frequency required to match the experimental 


results (Sec. 3.3.2). Thus, electron-neutral collisions are neglected. Because v ee and vn do 
not contribute to the momentum exchange between electron and ion fluids. [T2] collision 
frequencies for electrons and ions species are v e = v e i = 2.91 x 10 -12 n e j. e 


12 '" T e 3// "lnA and 




Vi = Vi e = m.eTOj V e j, respectively, from which we can see Vi e <C v e i. Here, T e and n e 
are given in eV and m~ 3 , respectively, and the Coulomb logarithm is calculated to be 
InA = 12. Singly ionised argon ions are assumed in this study, so that qi = —q e = |e|. 

The static magnetic field is assumed to be axisymmetric with B$ r <C Bq z and Bqq = 0. 
It is therefore appropriate to use a near axis expansion for the held, namely B$ z is only 
dependent on z and 

B 0 r {r, z) = ( 3 - 5 ) 


The antenna, as described in Sec. 1.2.2, is a left hand half-turn helical antenna. We assume 
that the antenna current is divergence free, to eliminate the capacitive coupling. Fourier 
components of the antenna current density are given by: 


jar — o, 


3 06 = Ia £lrn 2 1 S( r - Ra ) “ Za) + 5{z - Z a ~ L a )\ 

| H(z—Zq)H(z a +L a —z) imirW+tz—z„.)/L n \ \ 

L a ) 5 

_ j e —irmr[l + (z — z a )/La\ l_ e im -k , 

3a z = la 2 _ Ka > 

xH(z - z a )H(z a + La - z), 


(3.6) 


where L a is the antenna length, R a the antenna radius, z a the distance between the 
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antenna and the endplate in the source region, and H the Heaviside step function. Note 
that the antenna geometry selects only odd harmonic mode number m, as indicated by 
Chen et ah. m 

Equations B and ( |3.2[ ) are Fourier transformed with respect to the azimuthal angle 
and then solved (for an azimuthal mode number m) by a finite difference scheme on a 2D 
domain (r, z). The computational domain is shown in Fig. 3.2 The enclosing chamber 
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Figure 3.2: Computational domain employed to simulate the experimental setup shown in 
Fig. |1.8[ Here all dimensions are given in millimetres. 

is assumed to be ideally conducting so that the tangential components of E vanish at the 
chamber walls, i. e. 

Eg(T r ,z) — E z (L r , z) — 0, 

E r {r ,0) = Eg (r,0) = 0, (3.7) 

E r {r , L z ) = Eg(r, L z ) = 0, 

with L r and L z the radius and the length of the chamber, respectively. Moreover, all field 
components must be regular on axis, thus, Bg\ r= o = 0 and (rEg i)| r =o = 0 for m = 0; 
E z \ r=0 = 0 and (rEg )| r= o = 0 for m / O.fTO] In the present work, we choose the funda¬ 
mental odd mode number m = 1, which is preferentially excited in the helicon discharge 
launched by a left hand half-turn helical antenna. {7411751 129] . In the experiment, there is 
a radial air gap ( 0.055 < r < 0.0585 m) between the antenna and the glass tube, which is 
taken as glass region in the computational domain. We found that simulated results are 
insensitive to the dielectric constant in the glass region 0.05 < r < 0.055 m by varying the 
constant from 1 to 10. As no change was detected in the wave held, we therefore expanded 
the glass area radially to fill this air gap. The thickness of the antenna is approximately 
0.002 - 0.003 m. 


3.3.2 Computed and measured wave fields 

Based on the measured held strength configuration and plasma prohles shown in Fig. |3.1 


simulations are performed. Figure 3.3 shows the axial prohles of the computed B r ampli¬ 
tude and phase on axis, and their comparisons with experimental data. With the collision- 
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ality set to v e g = z/ e j, where v e $ is the effective collision frequency, and the antenna radius 
set to match the experiment, the predicted wave field is ~ 30% of the measured value, 
and the profile a poor match to the experiment. It is possible to obtain better agreement 
by varying the collisionality, which strongly affects the profile but leaves the magnitude 
largely unchanged, and the antenna radius, which strongly affects the field amplitude and 
leaves the radial and axial profiles of B unchanged. A qualitative match between mea¬ 
surement and simulation of the axial variation of B r is found using an enhancement in 
collisionality of = g e n(^ei + v ie) ~ Qen^ei with £ en = 9.5, and an adjustment in antenna 

dimension of i? s im = <ad^exp with t a d = 0.88. Since the collision frequency is proportional 
—3/2 

to T e , if the electron temperature was wrong, it has to be 4.5 times smaller in order 
to get the agreement. Calculation of the axial gradient of the computed phase variation 
shows a travelling wave, with a good agreement with data. 
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Figure 3.3: Variations of magnetic wave field in axial direction (on-axis): (a) |.B r | rms , (b) phase 
of B r . Computed results (lines: dotted for is e g = v ei and R slm = i? exp , dashed for z/ e g = v ei and 
-Rsim = 0.88I?exp, and solid for i/ e g = 9.5 v e i and R slm = 0.88i? eX p) are compared with experimental 
data (dots). 

The local minimum observed both experimentally and numerically in the axial profile of 
|B r | rms around z = 0.27 m has been also observed in many other devices. [81 )11301 1751173] for 
both uniform and non-uniform field cases. For the uniform field case, it has been suggested 
that the spatial modulation of the helicon wave amplitude is not caused by reflections from 
the end boundaries, but by a simultaneous excitation of two radial modes. @El EH ES [73j 
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Similarly, the minimum observed in MAGPIE cannot be explained by standing waves, 
because the amplitude becomes much smaller at bigger 2 (suggesting strong damping), 
and the phase advances with increasing z (denoting a travelling wave). Further, radial 


profiles of the wave field in Fig. 3.4 feature a possible superposition of the first and second 
radial modes of the m = 1 azimuthal mode. Therefore, we speculate that the minimum 
observed here may be also due to the simultaneous excitation of two fundamental radial 
modes. We will show later that radial gradient in plasma density is essential for the 
excitation of this local minimum under the present experimental conditions. 


(a) (6) 
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Figure 3.4: Variations of magnetic wave field in radial direction: (a), (c) and (e) are |B r | rms , 
l-B^Irms and |B 2 | rms , respectively; (b), (d) and (f) are the corresponding phase variations. Dots are 
experimental data while lines (dotted: z = 0.07 m, dashed: 0 = 0.17 m, solid: z = 0.21 m) are 
simulated results. 

Figure [3~4] shows the radial profiles of computed wave fields for v e R = 9.5^ and R s i m = 
0.88-Rexp at three axial positions in the target region, together with the experimental data 
measured at z = 0.17 m. The predicted wave field amplitude profile at z = 0.17 m is 
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consistent with the data, but the magnitude is nearly double the measured value, and the 
phase profile is a poor match. We have also computed the wave fields at axial locations 
with best agreement to the amplitude (z = 0.21 m) and phase (z = 0.07 nr). We justify 
this freedom of choice by the experimental uncertainty in axial density profile, and the 
numerical sensitivity identified in the radial profile of the wave field with axial position. 
Inspection of Fig. |3.4| reveals that it is possible to find a reasonable agreement to the 
wave amplitude and phase profile, albeit independently. As expected, all calculations 
show the wave mode structure of m = 1 through | B z (r = 0)| rms ~ 0, consistent with the 
antenna parity. The near-null minima in both |-Be(?’)| rms and |-B z (r)| rms suggest a likely 
simultaneous excitation of the first and second radial modes. This may account for the 
minimum in Fig. |3.3[ a), and is consistent with conclusions from others. [451 F73| 74, 75, 129 j 
Mori et al. m suggested that the superposition feature is associated with a wave focusing 
effect caused by the non-uniform magnetic held in the target region. 


3.4 Numerical profile scans 

It has previously been shown that the plasma density can be further increased by in¬ 
troducing a cusp or non-uniform static magnetic held in the vicinity of the helicon 
antenna. |18l mm f80lj To shed light on the increased plasma production, we perform 
a detailed numerical study on the effects of radial and axial plasma density gradients and 
axial magnetic held gradient on wave propagation characteristics. The enhancement of 
g en = 9.5 to u e i and the adjustment of t a d = 0.88 to R exp are still employed in this section 
because they provide a good agreement with the measured wave held. 


3.4.1 Axial profile of plasma density 


We hrst study the effect of varying the axial gradient of plasma density, which has been 
assumed to be linear with the static magnetic held so far, on wave propagations by com¬ 
paring the wave helds from three different on-axis density profiles shown in Fig. |3.5[ a). 
Other conditions are kept the same as previous sections. The computed wave helds in axial 


direction (on-axis) are shown in Fig. 3.5 b) and Fig. |3.5[ c). A log scale in the amplitude 
has been employed to see the wave propagation features clearly. We can see from the phase 
variations (Fig. |3.5[c)) that as density is decreased in the target region the wavelength 


increases, which is consistent with a simple theory developed previously (see Eq. (1.29)), 


3.83 


oj n e ep o 
k B n 


R p k Bq ^ 

Thus, if to, R p (plasma radius) and Bq are all hxed, k is proportional to n e , which means 
that the wavelength becomes larger at lower density. Here, the value of 3.83 is the hrst 
non-zero Bessel root of Ji(r) = 0, representing the hrst radial mode, which is assumed 
to be dominant in our case. Moreover, for all density profiles shown in Fig. |3.5[ a), the 
wavelength is bigger in region of 0 < z < 0.6 m than that in other regions, indicating a 
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locally increased phase velocity. Further inspection of Fig. |3.5| (b) and Fig. |3.5[ c) shows 
that density increasing in proportion to the static magnetic field has little effect on RF 
absorption, while the density level near the antenna affects the wave amplitude significantly 
at all axial locations. 


3.4.2 Axial profile of static magnetic field 


Second, following Sec. 3.4.1, we keep the axially uniform density profile and study the 
effects of axial gradient in static magnetic field, which is radially near uniform according 


to Eq. (3.5). Three employed field profiles are shown in Fig. |3.6[ a), which enable us 
to study the effects of field gradient in target and source regions seperately on wave 


propagations. Comparison between solid and dashed lines in Fig. 3.6 b) and Fig. |3.6|)c) 


shows that the axial gradient in magnetic field in the target region may increase the 
propagation distance of helicon waves, consistent with Mori et al.’s conclusion that a 
focused non-uniform magnetic field provides easier access for helicon wave propagations 
than a uniform held, m or perhaps the increased held strength in that region increases 


the collisional damping length of helicon waves. The simple theory in Eq. (3.8) is satished 


again here: with decreased held strength in the target region, the wavelength becomes 
shorter. Although the difference between dashed and dotted held prohles is small as shown 
in Fig. |3.6| (a) , the computed held amplitudes are signihcantly different. Figure [T6^ b) and 
Fig. |3.6[ c) show that with a uniform held prohle, the wave amplitude is much bigger 
than that with non-uniform held prohle for z < 0 m. Furthermore, waves keep their 
travelling features till the left endplate for uniform Bq, whereas for non-uniform Bq(z) the 
wavelength becomes smaller when approaching left, and the waves are not travelling at all 
when Bq(z) is low enough (z < —0.9 nr). 


3.4.3 Radial profile of plasma density 

Now, we keep the plasma density and static magnetic held both uniform in the axial 
direction, and study the effects of radial gradient in plasma density. The two density 
prohles employed are shown in Fig. |3.7[ a), with and without radial gradient, and their 
corresponding results are shown in Fig. E3 b ) and Fig. |3.7[ c). We can see first that 
the local minimum in wave amplitude prohles, e. g. at z = —0.52 m and z = 0.27 nr, 
disappear when the radial density prohle is hat, suggesting that the radial gradient in 
plasma density is essential to have a local minimum under the present conditions. Second, 
the wave amplitude is much bigger in both target and source regions for plasma density 
with radial gradient, suggesting that a radial gradient in density may be useful to nraxinrise 
the plasma production. 
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3.5 Collisionality and field direction 


3.5.1 Enhancement of electron-ion collision frequency 

In a similar manner to other work, uni ins] we have used an enhancement to v e i (here 
z/ e ff ~ 9.5is e i), in order to Hnd a qualitative match of simulated wave held to the data. In 
this section, we explore the physical consequences of scaling r' e f[/ iy ei hi simulations while 
keeping the adjustment t a( j = 0.88 in the antenna radius. Variations of wave amplitude 


on axis in the axial direction for different collision frequencies are shown in Fig. 3.8 As 
the collision frequency is increased from u e i to 15z/ e i, the wave amplitude decreases nearly 
everywhere, and the wave decay length is shortened. This indicates that the wave energy 
or power coupled from the antenna to the core plasma drops as the collision frequency 
becomes higher, and the power is more absorbed under the antenna. This is consistent 
with a previous conclusion that the RF energy is almost all absorbed in the near region of 
the antenna rather than in the far region. m The oscillations near the downstream and 
upstream ends at low i/ e g are caused by reflections from the ideally conducting endplates, 
which disappear if the endplates are moved further away. 

As suggested by Lee et al., [53] an enhanced electron-ion collision frequency may be 
due to ion-acoustic turbulence which can happen if the electron drift velocity exceeds the 
speed of sound in magnetised plasmas. Based on the experimental conditions in MAGPIE, 
we have calculated the threshold held strength B t, below which ion-acoustic turbulence 
can happen. This threshold is given by vd > C Sl where vd ~ kBT e /\e\BoR p is the electron 
drift velocity with ks Boltzmann’s constant and C s = yjksTe/rrii the speed of sound in 
magnetised plasmas. Application to MAGPIE conditions yields B x < 0.0224 T. Thus, the 
whole source region which produces helicon plasmas and waves, has a magnetic held below 
this threshold. The ion-acoustic turbulence has the effect of providing additional electron- 
ion collisions within a dielectric tensor model, and thereby improves the agreement with 
observations. Other possible reason for this enhanced electron-ion collisionality includes 
kinetic effects which are beyond the reach of the employed cold-plasma model, e.g. Landau 
damping. The Trivelpiece-Gould mode, which could be present in the low held region 
(B 0 (z) < 0.01 T) and is consistent with the strong edge heating observed at enhanced 
collision frequencies. [1311 132] is included in the present full wave simulations. For the 
present conditions of MAGPIE, the electron drift velocity is well blow the electron thermal 
speed, therefore, the Buneman instability does not occur. 


3.5.2 Direction of static magnetic field 

Observations have been made previously that the directionality of helicon wave propa¬ 
gations is dependent on the direction of static magnetic held in helicon discharges using 
helical antennas, but all in uniform held configurations. [73. 45,1, 1133 . [53] In this section, we 
study the directionality in a non-uniform held conhguration for a left-hand half turn heli¬ 
cal antenna. Specifically, we have computed the wave amplitude and wave energy density 
in MAGPIE for the experimental and held reversed configurations. In MAGPIE, the held 
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points from target to source. Figure 3.9 shows the computed axial profiles of wave ampli¬ 
tude on axis and 2D contour plots of wave energy density for both field direction pointing 
from target to source (Fig. |3.9| a) and Fig. |3.9[ b)) and field direction pointing from source 
to target (Fig. |3.9| and Fig. |3.9[ d)). In this calculation, we have chosen v e g = v e i to see 
more details, and chosen the density profile to be linear with Bq(z) in the axial direction 
and non-uniform in radius as measured in Fig. |3.1[ b). The field strength profile used here 
is shown in Fig. |3.1[ a). Inspection of Fig. |3.9| reveals that the wave energy is larger on the 
opposite side of the antenna, relative to the direction of the static magnetic field. This 
observation has been confirmed experimentally through finding that the plasma is brighter 
on the opposite side of the antenna relative to the direction of the applied external field. 
In summary, the dependence of the direction of helicon wave propagations to that of static 
magnetic field still exists even when the field configuration is non-uniform. 


3.6 Conclusions 

A RF field solver based on Maxwell’s equations and a cold-plasma dielectric tensor is 
employed to describe the wave phenomena observed in a cylindrical non-uniform helicon 
discharge, MAGPIE. Here, the non-uniformity is both radial and axial: the plasma density 
is dependent on r and z, the static magnetic field varies with z, and the electron tempera¬ 
ture is a function of r. A linear fitting was conducted for radial profiles of plasma density 
and electron temperature, and the fitted profiles were utilised in wave field calculations. 
A linear relationship between the axial profile of plasma density and the static magnetic 
field was also assumed. Other conditions used in the simulation were taken from experi¬ 
ment directly, including filling gas (argon), antenna current of 38.8 A, driving frequency 
of 13.56 MHz, and a left hand half-turn helical antenna. 

With an enhancement factor of 9.5 to the electron-ion Coulomb collision frequency i / e j 
to approximate the observed attenuation, and a 12% reduction in the antenna radius to 
match the amplitude of the wave field, the wave solver produced consistent wave fields 
compared to experimental data, including the axial and radial profiles of wave amplitude 
and phase. A local minimum in the axial profiles of wave amplitude was observed both 
experimentally and numerically, agreeing with previous studies. [8T1 11301 1751 173] Mode 
structure of m = 1 is consistent with the left hand half-turn helical antenna being used. 
A possible explanation for the enhanced electron-ion collision frequency has been offered 
through ion-acoustic turbulence, which can happen if the electron drift velocity exceeds 
the speed of sound in magnetised plasmas. [53] By calculating these two speeds based on 
MAGPIE conditions, we found that it is indeed satisfied in the source region of MAGPIE 
where the helicon plasmas and waves are produced. For the field strength under the 
driving antenna, the ratio of electron drift velocity and the speed of sound is around 2. 
Other possible candidate explanation may also include kinetic effects which are neglected 
in the cold-plasma model employed here. 

A numerical study on the effects of axial gradients in plasma density and static mag- 
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netic field on wave propagations was carried out. This showed that density increasing in 
proportion to the static magnetic field has little effect on RF absorption, while the den¬ 
sity level near the antenna affects the wave amplitude significantly at all axial locations. 
The axial gradient in magnetic field increases the decay length of helicon waves in the 
target region. The relationship between plasma density, static magnetic field and axial 
wavelength is consistent with a simple theory developed previously. [181] 

A numerical scan of the enhancement factor to u e i reveals that with increased electron- 
ion collision frequency the wave amplitude is lowered and more focused near the antenna. 
This is mainly because of stronger edge heating at higher collision frequencies which pre¬ 
vent more energy transported from the antenna into the core plasma. The wave amplitude 
profile at v e R = 9.5 v e i, which agrees with experimental data, shows consistent feature with 
a previous study that the RF energy is almost all absorbed in the near region of the an¬ 
tenna rather than in the far region. m We also studied the effect of the direction of 
static magnetic field on wave propagations, and found the antiparallel feature that waves 
propagate in the opposite direction of magnetic field for the antenna helicity existing in 
these experiments. This dependence of the direction of helicon wave propagations to that 
of static magnetic field in a non-uniform field configuration is consistent with previous 
observations made in uniform field configurations. [711 S5 111331 [53] 

Physics questions raised by this work include: further explanation of exactly how 
axially non-uniform field might affect the radially localised helicon mode, [35] inclusion of 
different m numbers in the glass layer and any subsequent coupling to the plasma at the 
plasma-glass interface, and identification of independent first and second radial modes that 
superpose to yield a local minimum in wave field amplitude at z = 0.27 m. Experimental 
measurements that might corroborate the wave field generation mechanism and associated 
physics include: the measurement of axial profile of density, and measurements of n e , B, 
and E with a reversed field. 
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Wave propagations in a pinched plasma 


(a) 



z (m) 


Figure 3.5: Normalised axial profiles of plasma density and the corresponding variations of mag¬ 
netic wave field in axial direction: (a), normalised axial profiles of plasma density, (b) log scale of 
I B r |rms, and (c) phase of B r . They are all measured on axis. In (a), the solid line is associated 
with density profile linear with Bq(z) : the dashed line is the same to the solid one except in region 
of 0 < z < 0.7 m where the density maximum is adjusted to match preliminary experimental 
observations, and the dotted line represents a z-independent density profile. 
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Figure 3.6: Axial profiles of static magnetic field and the corresponding variations of magnetic 
wave field in axial direction: (a) axial profiles of static magnetic field, (b) log scale of |£? r | rms 
(on-axis), and (c) phase of B r (on-axis). In (a), the solid line shows original experimental data 
(Fig. |3.1| (a)), while the dashed line shows the same except being flattened in region of z > 0 m 
and the dotted line is flattened everywhere. 
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Wave propagations in a pinched plasma 


(o) 



r (m) 


(b) 



(c) 



Figure 3.7: Radial profiles of plasma density and the corresponding variations of magnetic wave 
field in axial direction: (a) radial profiles of plasma density, (b) log scale of |.B r | rms , (c) phase of 
B r . All parameters are measured on axis. 
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Figure 3.8: Variations of on-axis wave amplitudes in axial direction with different electron-ion 
collision frequencies: log scale of B r | rms • 
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Wave propagations in a pinched plasma 
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Figure 3.9: Axial profiles of magnetic wave field (on-axis) and contour plots of wave energy 
density in (z, r) space for a non-uniform plasma density: (a), log scale of |-B r | r ms for upstream 
B 0 {z)\ (b), wave energy density for upstream B 0 (z ); (c), log scale of |i?r.| r ms for downstream B 0 (z); 
(d) wave energy density for downstream Bq(z). 






































Chapter 4 


Gap eigenmode of radially 
localised helicon waves in a 
periodic structure 


The EMS (ElectroMagnetic Solver [123] ) is employed to model a spectral gap and a gap 
eigenmode in a periodic structure in the whistler frequency range. A RLH (Radially Lo¬ 
calised Helicon(35]) mode is considered. We demonstrate that the computed gap frequency 
and gap width agree well with a theoretical analysis, and find a discrete eigenmode inside 
the gap by introducing a defect to the system’s periodicity. The axial wavelength of the 
gap eigenmode is close to twice the system’s periodicity, which is consistent with Bragg’s 
law. Such an eigenmode could be excited by energetic electrons, similar to the excita¬ 
tion of Toroidal Alfven Eigenmodes (TAE) by energetic ions in tokamaks. Experimental 
identification of this mode is conceivable on the LAPD.[S] 


4.1 Introduction 

It is a generic phenomenon that spectral gaps are formed when waves propagate in peri¬ 
odic media, and eigenmodes can exist with frequencies inside the spectral gaps if a defect 
is introduced to break the system’s perfect translational symmetry, and create an effective 
potential well to localise these waves. [T34 . 1135 . 1136111371113811139] Fusion plasmas have 
a few periodicities that can produce spectral gaps: geodesic curvature of held lines. (140 ] 
elongation [TJX, 142 ] or triangularity of flux surfaces .( 113) helicitv [l44l and periodic mirror¬ 
ing in stellerators. [H5j Weakly damped eigenmodes which are readily destabilised by ener¬ 
getic ions often reside in these gaps, and they may degrade fast ion confinement. [146 . 147 ] 
The most extensively studied gap eigenmode is TAE. (148] however, there are also numer¬ 
ous other modes with similar features. (1491 [150 1 11511 [152] 

Zhang et. al. [10] observed a spectral gap in the shear Alfven wave continuum in experi¬ 
ments on the LAPD with a multiple magnetic mirror array, and obtained consistent results 
through two-dimensional numerical modelling using the finite difference code, EMS. (124 ] 
Although eigenmodes inside this gap were not formed, a possible experimental implemen¬ 
tation was proposed to detect them. The idea is to use the endplate of the machine as a 
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Gap eigenmode of radially localised helicon waves in a periodic structure 


“defect” that breaks the axial periodicity or to vary current in one of the independently 
powered magnetic coils. 

In this chapter, we will use EMS to examine a gap eigenmode in a linear system with 
slightly broken axial periodicity. We will consider the RLH mode, whose radial structure 
has been described in m- Beyond [35], which deals with uniform magnetic field, the 
spectral gap and gap eigenmode of the RLH mode will be theoretically and numerically 
studied in a rippled magnetic field in the present study, together with those for shear Alfven 
waves in Chapter [5j We will show that the computed gap frequency and gap width agree 
well with the theoretical analysis, and that there is a discrete eigenmode inside the gap. 
Such an eigenmode could be excited by energetic electrons, similarly to the excitation of 
TAE by energetic ions in tokamaks, and possibly observed on a linear device, e. g. LAPD. 


4.2 Theoretical analysis 

4.2.1 Basic equations 


The spatial structure of RLH wave field and the corresponding dispersion relation can 


be found from Eq. (1.52) that applies to a whistler-type linear wave in a cold plasma 
cylinder. m For a cylindrical coordinate system of (r, (p , z) with ip the azimuthal angle 
and wave form of exp[z(m<^ + k z z — ut)], it has the form 


1 d dE m 2 m u 2 EdD/dr 
r dr dr r 2 k 2 r c 2 1 — {mdD/dr)/k 2 rP 


where E = E z — (k z r/m)E v with E z and E p the axial and azimuthal components of 
the wave electric field, respectively, D = —u> 2 /lolo c , and P = —uj 2 /uj 2 . Here, oj p and u c 
are the electron plasma frequency and electron cyclotron frequency, respectively. In this 
section, we will limit our consideration to the case of sufficiently dense plasma in which 
uj 2 3> c 2 m/a 2 , where a is the plasma radius, and nonzero values of the azimuthal mode 
number m as the radial nonuniformity of plasma density has a surprisingly strong effect 
on the structure of helicon modes with m ^ 0. [35] The resulting eigenfrequency for the 
mode of interest scales roughly as[35] 


u 


u, 


2 c 2 

2 ' 


kic 


u, 


(4.2) 


which shows that the term — {mdD/dr)/k 2 rP in the denominator of Eq. (4.1) can be 
estimated as 


— {mdD/dr)/k 2 z rP ~ m !/ , C 1. (4.3) 

a u p 

It was also shown in [35 j (see Eq. (16) there) that 

E z —{mdD/dr)/k 2 rP 
E 1 — (mdD/dr) jk 2 rP 


(4.4) 
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which is thus much less than unity for the mode of interest. The underlying reason is 
that the electron conductivity along the guiding magnetic field is much greater than the 
cross-field Hall conductivity. It is therefore allowable to set E~ = 0 and ignore the term 


— (mdD / dr) / k 2 rP in the denominator of Eq. (4.1), which gives 

7 v i d d 2n \ u 2 dD 

fc, —r—rE, n — m E m = E^m—^r 


dr dr v v ^ c 2 dr 


(4.5) 


We now generalise this equation to the case of slightly modulated (^-dependent) plasma 
equilibrium with a following separable form of cz 2 /u) c 


P w po( r )ri f \ ( M 

— =-[1 — e(z) cos(g, 2 )J. 

w c o 


(4.6) 


Here, t < 1 and q are the modulation envelope and wavenumber respectively. The axial 
scale-length of the envelope e(z) is assumed to be much greater than 1/q. We recall that 


k 2 = —d 2 /dz 2 and transform Eq. (4.5) to 


dz 2 


{dr r dr rEv ~ m2Ev ) = Eq>m ^ [1 ~ C0 ^ ( l z )\ r l r • (4-7) 


The modulated equilibrium introduces resonant coupling between the modes with k z = q/2 
and k z = —q/2, which suggests the following form for E v 


Ep = A + e iqz ! 2 + A_e~ iqz ' 2 , 


(4.8) 


where A + and are slow functions of z compared to cos(gz). Let ujq be an eigenfrequency 


of Eq. (4.7) for e = 0 and \k z \ = q/2. It is then straightforward to separate spatial scales 


in Eq. (4.7) and obtain a set of coupled equations for A + and 


(wv r iv rA + ~ m2A +) - A + m! f r w (So) 


- ni?A + ) + A +m ^rf r (§) - 


(4.9) 


( W r Wr rA ~ ~ m2A ~) ~ A - m ^ r §r (So) = 


+*4 {lrlrA--m 2 A.)+A.m^rl (g) - \A + m^e{z)rl (g) 


(4.10) 


We have neglected second axial derivatives of A + and and arranged Eq. (4.9) and 


Eq. (4.10) so that their right hand sides represent relatively small terms compared to the 


dominant terms on the left hand sides. With this ordering, we conclude that the radial 
dependencies of A + and A _ need to be close to the eigenfunction T(r) of the lowest order 
ODE 
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Gap eigenmode of radially localised helicon waves in a periodic structure 


The differential operator in this equation is self-adjoint. As a result, multiplication of 


Eq. (4.9) and Eq. (4.10) by rT and integration over radius lead to: 


0 = ~ i( lWz f rq> (m r m rA + ~ m 2 A + ) dr + / rfd+m^/rl dr 


~h f rVA-m^e(z)r£ dr, 


dr yu c0 

0 = i( i§z / r ^ ( lk r lk rA ~ ~ m 2A -) dr + J r'fyA-m^^-rdr 

-\jrA>A + m^e(z)r-§- dr. 

Having eliminated the lowest order terms, we now set: 

A + = F(z)A?(r), A- = G(z)A/(r) 


(4.12) 


(4.13) 


(4.14) 


in Eq. (4.12) and Eq. (4.13) to obtain the following set of coupled equations for F(z) and 
G(z): 

Ai f)T? (,) — Mr* Ct 

(4.15) 


4 i dF cjj — ojq G 

0= + F- 

q oz ujq 2 


AidG u-uj 0 F 

0 =- — +G - yz£\z)- 

q oz wo 2 


(4.16) 


4.2.2 Spectral gap and continuum 


In the case of ^-independent e, Eq. (4.15) and Eq. (4.16) have exponential solutions 


F oc G oc e* 


(4.17) 


with k the wave number. The corresponding dispersion relation has two roots: 

u + ( k ) = W 0 [l + v / ( 4k /<?) 2 + (e/2) 2 ], 

oj~{k) = w 0 [l - \/(4 n/q) 2 + (e/2) 2 ]. 


(4.18) 


Here, w+(k) and w_(k) are the continuum frequencies above and below the spectral gap, 
respectively. The upper and lower tips of the spectral gap correspond to k = 0: 


w+(0) — ujo (l + |) , 
u;_(0) = ujq (1 - §) . 


(4.19) 


The normalised width of the gap, Aw = [w+(0) — w_(0)]/wo, is therefore equal to the 
modulation amplitude e. The central frequency of the gap wq needs to be found from 


Eq. (4.11) as the RLH eigenfrequency for an axially uniform plasma cylinder and it can 
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be written as 


w 0 


w c0 (0)cV 
4w p0 (0) 2 ’ 


(4.20) 


where w c o(0) and w p o(0) are the on-axis values and the numerical form-factor T is deter¬ 
mined by the plasma radial profile. 


4.2.3 Wall-localised eigenmodes 

A discrete-spectrum eigenmode can be created inside the spectral gap if the system’s 
periodicity is broken. We illustrate that by considering an ideally conducting endplate 
located at z = zo, so that the plasma now occupies only a half cylinder to the right of the 
endplate (z > zq). The boundary condition at the endplate is 


E<p(r, zo) = 0 


(4.21) 


or, equivalently, 


F(z 0 )e iqzo/2 + G{z 0 )e~ iqzo/2 = 0 


(4.22) 


(see Eq. (4.8) and Eq. ( 4. 14| ) ). The electric field of the discrete-spectrum mode must also 
vanish at z —> oo, i. e., 

F( oo) = G(oo) = 0. (4.23) 


Equation (4.15) and Eq. (4.16) (with e = const ) admit an exponential solution 

F oc G oc e~ Xz 


(4.24) 


with a positive value of A that satisfies boundary conditions in Eq. (4.22) and Eq. (4.23). 


Indeed, the exponential ansatz, Eq. (4.24), reduces Eq. (4.15), Eq. (4.16) and Eq. (4.22) 
to an algebraic set: 


4iA _ ^w — wn G 

0 =- F + F -- — —e, 

q w 0 2 

o = ia a + g u-^_F 
q w 0 2 


F + Ge~ iqz ° = 0. 

The solvability conditions for this set determine w and A as functions of zq as follows: 


(4.25) 

(4.26) 

(4.27) 


-—— = ~\ cos(qz 0 ), (4.28) 

Wo 2 

A = lqsin(qz () ). (4.29) 

O 

We observe that A is positive for 0 < qzo < vr, and the wave held therefore vanishes at 
infinity in this case to represent a discrete eigenmode whose frequency is automatically 
inside the gap. In particular, the eigenfrequency is exactly at the gap centre when qzo = 
ir/2. We further note that this wall-localised solution can also be viewed as an odd-parity 
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eigenmode (E (fi (z) = —E v (2zq — z)) in the entire periodic cylinder with a defect at z = zo. 
An axial profile of Wp/w c with such defect is shown in Fig. 4.1 'a). 


(a) 



(b) 



Figure 4.1: Axial profiles of Wp/w c with defects: (a) defect location at sin(gzo) = 1, (b) defect 
location at sin(gzo) = — 1. 


4.2.4 Even-parity eigenmode 

In contrast with Fig. |4.1[ a), Fig. |4.1[ b) illustrates a defect with ir < qzq < 2i r that does 
not produce an odd-parity mode. Instead, an even-parity mode exists in this case. This 


mode is still described by Eq. (4.15) and Eq. (4.16), but the boundary condition at z = zq 


now changes from Eq. (4.21) to 


dE v (r, z) 
dz 


= 0 , 


J z=zo 


or, equivalently, 


F - Ge~ iqz ° = 0. 


The solvability conditions for Eq. (4.25), Eq. (4.26) and Eq. (4.31) now give: 


uj -uj 0 e 

— = -cos(gz 0 ), 

o; 0 2 


(4.30) 


(4.31) 


(4.32) 


A = --qsm(qzo), 


(4.33) 


and we observe that A is positive (the wave field vanishes at infinity) when sm(qzo) < 0, 
i. e. this even-parity eigenmode indeed requires n < qzo < 27r for its existence. 


4.2.5 Schrodinger equation for gap eigenmodes 


Equation (4.15) and Eq. (4.16) can be straightforwardly transformed into two independent 


second-order equations for F — G and F + G: 

16 8 1 d(F-G) 


q 2 dz ( u-uq _ dz 

( «0 2 J 


+ (f _ g)( ^2 + 1 ] = °, 


(4.34) 
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16 d 


q 2 dz ( u-WQ I dz 

V WO + 2) 


a( - F+a '> + (F+G) lu - uo 


w 0 


= 0. 


(4.35) 


Both of them can be further reduced to a time-independent Schrodinger equation when 
e is nearly constant. We assume e = (e) + u(z) with u(z) <C (e) and u(±oo) = 0. The 
discrete-spectrum modes in this case are very close to the tips of the gap. For the lower 
tip, we have 

t, ) - /, In 11=\ 

(4.36) 


w-w 0 (e) 

- = ~T + S - 


w 0 


with |d_| <C (e). We can then neglect u(z) and 5- in the derivative term of Eq. (4.34) and 
get 




dz 2 ' 16 

Similar procedure for the upper tip gives 

d 2 (F + G) , q 2 (e) 


8- + 


u(z) 


= 0. 


dz 2 


+ hr (F+G) 


-<*+ + 


u{z) 


= 0. 


(4.37) 


(4.38) 


We note that any negative function u(z) acts as a potential well that supports spatially 


localised eigenmodes of Eq. (4.37) and Eq. (4.38) with 5- > 0 and 5+ < 0, respectively. 


The eigenfrequencies of these modes belong to the spectral gap in the continuum. 


4.3 Numerical results and discussion 


The EMS. [124j as introduced in Sec. 3.3.1 is employed to study the RLH spectral gap and 


gap eigenmode inside. The computational domain is shown in Fig. 4.2 A helical antenna 
is employed to excite an m = 1 mode in the plasma. The enclosing chamber is assumed to 


be ideally conducting so that the boundary conditions are the same to those in Eq. (3.7). 


nr 
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Figure 4.2: Computational domain. The solid bar denotes a half-turn helical antenna. The dot- 
dashed line is the machine and coordinate system axis (r = 0). The coordinate system (r, ip, z) is 
right-handed with an azimuthal angle ip. 
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Gap eigenmode of radially localised helicon waves in a periodic structure 


4.3.1 RLH mode in a straight cylinder 

We first use EMS to calculate the radial structure of the RLH mode in a straight cylinder 
with a uniform static magnetic field and recover the mode dispersion relation. We con¬ 
sider a single-ionised argon plasma with a radial density profile shown in Fig. |4.3[a). The 



k, (m ') r (m) 


Figure 4.3: Plasma density profile and RLH wave field in a straight cylinder: (a) radial profile 
of unperturbed plasma density, (b) azimuthal magnetic field of the m = 1 mode for an illustrative 
frequency cu = 31.4 Mrad/s (solid bar shows antenna location, and the endplates are located at 
z = 0 m and z = 30 nr respectively), (c) computed dispersion relation for the m = 1 mode (dots) 
and analytical scaling from Eq. ( |4.20| ) with q/2 = k z and T = —2.03 (dashed curve), (d) radial 
profile of the m = 1 mode at z = 20 m for ui = 31.4 Mrad/s (the profile shows some modulations 
due to coupling to other weaker modes produced by the antenna). 


surrounding vacuum chamber and the position of the rf-antenna are shown in Fig. 4.2 


The static magnetic field in this simulation is Bq = 0.01 T and the electron-ion collision 
frequency on axis is z/ e j| r= o = 4 x 10 6 s -1 . We run EMS for various frequencies of the 
antenna current, select the m = 1 azimuthal component of the plasma response and cal¬ 
culate the dominant axial wavenumber, k~, of this component via Fourier decomposition. 
We observe an apparently preferred axial coupling direction for the m = 1 wave generated 
by the helical antenna. The collision frequency is sufficient to damp the radiated waves 
significantly towards the endplates but the dominant wavelength is still seen clearly in the 
plot of the excited field (Fig. |4.3| (b)). The resulting relation between k z and oj is shown 
in Fig. |4.3| (c). We note that k z scales as y/uj with frequency, in agreement with analytical 
scaling obtained in [35], namely ui ~ u: c k z c 2 /ujp. A least squares fit to the dispersion rela- 
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tion given by Eq. (4.20) with q/2 = k~ shows that the form factor T equals —2.03 for the 
selected density profile. The radial profile of this m = 1 mode is presented in Fig. |4.3[d). 


4.3.2 Spectral gap in a periodic structure 

In order to demonstrate the spectral gap, we run EMS for the same conditions as in 


Sec. |4.3.1| except that the static magnetic field is weakly modulated along z and has the 
form Bq z = Bq[1 + 0.25cos(g,?)] with Bq = 0.01 T and q = 40 m _1 . The radial component 


of the static field is determined by Eq. (3.5). We scan the antenna frequency (for a fixed 


amplitude of the antenna current) and analyse the wave propagation from the antenna. 


Figure 4.4 illustrates the results. Figure |4.4[ a) is a plot of the on-axis signal at a 4.8 m 
distance from the antenna, and shows an interval of strong suppression (from 36 Mr ad/s to 
44 Mrad/s), which represents the anticipated spectral gap. The location and width of this 
gap agree well with analytic estimates based on Eq. (4.19) and Eq. ( 4.20[ ). Figure |4~4) (b) 
presents spatial distributions of wave energy for frequencies that are inside and outside of 
the spectral gap, from which we also see that wave propagation is evanescent inside the 
gap. 


4.3.3 Gap eigenmodes 

To form an eigenmode inside the spectral gap shown in Fig. |4.4[ we introduce a defect 


in the otherwise periodic static magnetic field. Equation (4.28) and Eq. (4.32) suggest 


that this defect should be located at cos(qzo) = 0, in order for the eigenmode frequency 
to be at the gap centre. The mode is then expected to have the shortest possible width 


(see Eq. (4.29) and Eq. (4.33)). Moreover, Sec. 4.2.3 and Sec. 4.2.4 indicate that the gap 


eigenmode can have either odd parity or even parity, depending on the defect profile. We 
will confirm this analysis by simulation results shown in Fig. |4.5| and Fig. |4.6| 

Figure |4~5| (a)-(d) illustrates the odd-parity gap eigenmode associated with the defect 
shown in Fig. |4.1[ a). In order to minimise the role of collisional dissipation and thereby 
obtain a sharper resonant peak in Fig. |4.5[c), we have reduced the collision frequency to 


0.1 v e i in Fig. 4.5 simulation results. Figure |4~5)(a) shows that the eigenmode is a standing 


wave localised around the defect. Its exponential envelope with a decay length of 1.25 m is 


consistent with analytical expectation from Eq. (4.29). The evanescent held of the mode is 


only weakly coupled to the distant antenna, but the resonance with the antenna frequency 
still allows the mode to be excited easily. The internal scale of this mode is close to twice 


the system’s periodicity, consistent with Bragg’s law. Figure 4.5 b) shows the on-axis 
profile of the mode electric held, implying an odd function of E, p at zq. The corresponding 
magnetic held B^ is, instead, an even function. Figure [475| (d) presents a full view of the 
plasma response in (z, ui) space with a clear eigenmode peak inside the spectral gap. 

We have also performed similar calculations for the defect shown in Fig. |4.1| (b). This 
defect produces an even-parity mode seen in Fig. |4.6[ a)-(d). Except for different symme¬ 
tries, the features of the even and odd modes are apparently similar. 
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0 5 10 15 20 

z(m) 


Figure 4.4: Identification of the spectral gap in simulations: (a) on-axis wave field at 2 = 15 m as 
a function of driving frequency; vertical lines show analytical predictions for the gap centre uiq and 
the upper and lower tips of the continuum and w_, (b) spatial distributions of wave energy for 
antenna frequencies inside and outside the spectral gap (solid bar denotes the antenna location). 


The effects of collisionality on gap eigenmodes are shown in Fig. 4.7, from which we 
could see clearly a strength drop as the collisionality is increased. Although the resonant 
peak decreases and broadens for larger values of u e i. it is still clearly visible even at the 
highest collision frequency. 


4.4 Summary and thoughts about possible experiments 

We have shown that longitudinal modulation of the guiding magnetic field in a plasma 
column creates a spectral gap for radially localised helicon waves. Calculations performed 
with an EMS code reveal that this gap prohibits wave propagation along the column when 
the driving frequency of the RF antenna is in the forbidden range. The calculated width 
of the gap is consistent with analytical estimates. We have also shown that a discrete 
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(a) 


(b) 
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Figure 4.5: Odd-parity gap eigenmode: (a) longitudinal profiles of the static magnetic field (solid 
line) and the RF magnetic field on axis for lo = 40.2 Mrad/s (dots), together with the theoretically 
calculated envelope (dotted line), (b) longitudinal profile of E v (on-axis) at w = 40.2 Mrad/s 
(vertical clot-dashed line marks the defect location, and solid horizontal bar marks the antenna), 
(c) resonance in the dependence of the on-axis amplitude of the RF magnetic field on driving 
frequency at the location of the defect, (d) 3D plot of the on-axis wave held strength as a function 
of z and w. 


eigenmode can been formed inside the spectral gap by introducing a local defect to the 
periodic structure. Both the theoretical analysis and simulations demonstrate two types 
of gap eigenmode in the “imperfect” system: odd-parity mode and even-parity mode, 
depending the type of the defect employed. The gap mode is localised around the defect 
and represents a standing wave rather than travelling wave. Its distinctive feature is a 
resonant peak in the plasma response to the antenna current. The gap eigenmode has two 
characteristic spatial scales: a short inner scale and a smooth envelope. The inner scale 
is nearly twice the system’s periodicity, which is characteristic for Bragg’s reflection; the 
envelope depends on the modulation amplitude and it scales roughly as the inverse width 
of the spectral gap. 

A plausible way to identify the gap mode in a linear device with multi-mirror con¬ 
figuration would be to use the end-plate of the machine as a controllable defect in the 
periodic system. This could make the mode observable with a modest number of mir¬ 
rors in the machine. The LAPD is an apparent candidate for such experiments, provided 
that dissipative processes in the plasma do not destroy the gap-mode resonance. We have 
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Figure 4.6: Even-parity gap eigenmode: (a) longitudinal profiles of the static magnetic field (solid 
line) and the RF magnetic field on axis for ui = 41.8 Mrad/s (dots), together with the theoretically 
calculated envelope (dotted line), (b) longitudinal profile of E v (on-axis) at w = 41.8 Mrad/s 
(vertical dot-dashed line marks the defect location, and solid horizontal bar marks the antenna), 
(c) resonance in the dependence of the on-axis amplitude of the RF magnetic field on driving 
frequency near the location of the defect, (d) 3D plot of the on-axis wave field strength as a 
function of z and w. 


estimated these dissipative processes for LAPD plasma conditions. Collisional damping 
rate of the m = 1 RLH mode can be roughly estimated as u e c 2 /(a 2 LJp),\ 32] where v e is 
the collision frequency for electrons (a sum of electron-ion and electron-neutral collisions). 
Collisional damping can therefore be controlled by varying the plasma radius. Our illus¬ 
trative simulation is for a = 0.05 nr, but the LAPD vacuum chamber allows to produce 
plasmas with considerably larger value of a (up to 0.5 m[[f0] 0]) and thereby reduce the 
damping significantly. It may be preferable to choose plasma parameters in such a way 
that the gap mode frequency is in the range of operational frequencies for standard RF 
generators (6 — 28 MHz|18|) through Eq. (4.20). The LAPD magnetic coils are about 
0.61 nr in radius (inner),[9] which limits the number of magnetic field periods within the 
machine length to about 20. This number seems to be large enough for the mode identi¬ 
fication, especially because the axial extension of the mode is expected to decrease when 
the modulation amplitude increases. Stronger modulation of the magnetic field broadens 
the gap and makes it easier to identify the mode in the presence of hnite damping. Thus, 
the amplitude of modulation is an apparent control knob for the conceivable experiment. 
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Figure 4.7: Effects of collisionality on gap eigenmodes: (a) results for odd-parity mode 
(Fig. |4.5[c)), (b) results for even-parity mode (Fig. |4.6[c)). 


The LAPD experiment can use larger plasma radius and stronger magnetic field (with 
the same plasma density n e ~ 10 18 m~ 3 ) to reduce the damping but still keep the mode 
frequency in the convenient operational range for the RF-antenna. 
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Chapter 5 


Gap eigenmode of shear Alfven 
waves in a periodic structure 


In this chapter, a gap in the shear Alfven wave spectrum is predicted for a cold cylindrical 
plasma by modulating the magnetic field periodically. A discrete eigenmode is generated 
inside the gap by breaking the field periodicity through an introduced defect. We find that 
the axial wavelength of the gap eigenmode is nearly twice the system’s periodicity, which 
is characteristic of Bragg’s reflection, and the decay length of the eigenmode agrees well 
with the analytical estimate. Our treatment is based on numerical and analytical solutions 
to the cold plasma dielectric tensor, and complements Chapter [4] on the gap eigenmode 
generation in the whistler frequency range. 


5.1 Introduction 

In Chapter [4j a discrete eigenmode of RLH waves was generated inside a spectral gap by 
introducing a defect to the system’s periodicity. [15] This gap eigenmode is in the whistler 
frequency range and could be excited by energetic electrons, similar to the excitation 
of toroidal Alfven eigenmodes by energetic ions in tokamaks. We concluded that this 
eigenmode is observable on the LAPD, given that dissipative processes in the plasma do not 
destroy the gap-mode resonance. The LAPD, which has 10 sets of pancake electromagnets 
(90 in total) surrounding a 17.56-nr-long, 1-m-dianr plasma column, is capable of forming 
multiple magnetic mirrors and identifying a spectral gap of plasma waves. m This makes 
it a promising candidate to identify a gap eigenmode. However, to observe the RLH gap 
eigenmode, plasma parameters on the LAPD have to be chosen carefully to keep the mode 
frequency within the convenient operational range of standard RF generators. In 2008, 
Zhang et al. m observed spectral gaps and continua of SAW on the LAPD, with the gap 
frequency about 0.3 of the ion cyclotron frequency. Although eigenmodes were not formed 
inside this gap, they proposed a possible experimental implementation to detect them: 
varying current in one of the independently powered magnetic coils to introduce a strong 
defect that breaks the system’s periodicity. Therefore, we are interested in generalising the 
gap-mode analysis from the whistler frequency range to ion cyclotron range of frequencies. 
The SAW, which have frequencies below the ion cyclotron frequency and are fundamentally 
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different from those with frequencies above, will be focused on with reference to Zhang et 
al.’s experiment. [10] We shall numerically form a gap eigenmode of SAW by introducing 
a defect to the system’s periodicity. 


5.2 Theoretical analysis 

Although the focus of the study is on SAW, the theoretical analysis presented below is 
also valid for waves with frequencies slightly above the ion cyclotron frequency. 


5.2.1 Basic equations 


We consider an electromagnetic wave in a cold plasma cylinder with uniform static mag¬ 
netic held and radially non-uniform plasma density. The cold plasma approximation im¬ 
plies a much greater longitudinal wave phase velocity than particle thermal velocity, and 


assures low Landau damping. We start from the linear wave equation (see Eq. (1.44)) 


V x V X E = %-e ■ E. 


(5.1) 


For in « uj C i, the components of e become 


5 = 1- 


,2 

v pe 


it, 


i \ _ ^ce^p e ^ 

~ “(uj 2 -uj2 e ) T 


P= 1 - 




(5.2) 


Here, the relation Up e /u ; ce = —uj^/ud has been used. As in Chapter |dj we will limit our 
consideration to the case of sufficiently dense plasma and nonzero values of the azimuthal 
mode number m, because the radial nonuniformity of plasma density has a surprisingly 
strong effect on the structure of m / 0 modes. [35] Due to w < co pe , the parallel conduc¬ 
tivity P is very large so that the axial component of the electric held vanishes, namely, 


E z = 0. This simplifies the radial and azimuthal components of Eq. (5.1) as: 
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It is convenient to combine Eq. (5.3) and Eq. (5.4) into 
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and use Eq. (5.3) and Eq. (5.5) instead of Eq. (5.3) and Eq. (5.4). It follows from Eq. (5.3) 
that 

1 F) im 

0 (5.6) 


1 d , _ . im _ 
~ — (rE v ) - E r 


for modes with sufficiently small values of k z and u>. We choose small values of k z to see 
clear effects of the radial non-uniformity in plasma density on the wave mode structure 


for non-zero azimuthal mode numbers. [35] We can then use Eq. (5.6) to eliminate E r from 


Eq. (5.5), and obtain the wave equation of SAW 


0=[^S-k 2 z \rE v - 


rE, n u; 2 dD 


m c 


dr 


d 

+ r W 

dr 




(5.7) 


Equation (5.7) is identical to the wave equation of RLH waves when S' = 0 (see Eq. (1.52)) 


and can be solved in a similar way. 


5.2.2 Uniform magnetic field: dispersion relation 


We first solve Eq. (5.7) for a uniform magnetic field to derive the dispersion relation of 


SAW. We shall construct a surface-wave solution through considering a step-like radial 


profile of plasma density (see Eq. (1.54)) 


m = 


n_, r < ro 
n+, r > r 0 


(5.8) 


with ro the radius of the density discontinuity. The width of the discontinuity layer at 
r = ro is assumed to be larger than the skin depth, namely 2S W > c/ui pe with 5 W half 
width of the layer. This density profile leads to: 


W = t 5 ( re +) “ S(n_)]<5(r - r 0 ), 
fy = [D{n+) - D(n-)]5(r - r 0 ). 


(5.9) 


We then assume a specific form of E v (see Eq. (1.55)) 

(r/r 0 )H -1 , r < r 0 
(r/ro)“l m l~ 1 , r > r 0 


ro 


(5.10) 


with Eq a constant electric field. Equation (5.7) is satisfied at r / ro for any u and k 


for the density profile constructed. For modes caused by the density jump, we integrate 
Eq. (5.7) across the discontinuity layer ro — S w —> ro + 5 W 

(w* S ~ X ) r ^ rE F) . - Iro^t ~ X ) rr^Etpdr 


- r 0 +(5 w 
1 r 0 -<5 w 

= T&mrE v <gdr. 


(5.11) 
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The right hand side of Eq. (5.11) becomes 

,2 


UJ 

k 2 c 2 


mEo[D(n +) — D(n-)] 


(5.12) 


due to Eq. (5.9), and the second term on the left hand side vanishes because of the small 


value of |5 W | and finite Eq, 


rro+S^ 


/ , l U2 r 2 

'ro-dy, \^z L 


S — 1 | m 2 E l ndr = 0. 


(5.13) 


Substituting Eq. (5.2), Eq. (5.9) and Eq. (5.10), Eq. (5.11) can then be rewritten as 
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(5.14) 


Here, T = (w 2 j_ — uj 2 i+ ) / (u 2 ^ + uj 2 i+ ) with u p i- and oj p i + the ion plasma frequencies 
inside and outside of the discontinuity layer, respectively, and x = oj/u C i is the normalised 


wave frequency. Equation (5.14) describes the dispersion relation of SAW propagating 


in a plasma cylinder with uniform magnetic field and step-like radial profile of plasma 


density. For SAW with u <C oj c i and propagating in a uniform plasma, Eq. (5.14) becomes 


ua = fczi'A due to x < 1 and \k = 0. This is the well-known dispersion relation for SAW 
in a slab geometry. 


5.2.3 Periodic magnetic field: spectral gap and continuum 


We now introduce a slight periodic modulation to the static magnetic field, and solve 
Eq. (5.7) for the corresponding SAW spectrum. The modulated field is of the form B z = 
Bo [l + ecos((/, 2 )] with Bq the equilibrium field strength, e<l and q the modulation factor 
and wavenumber, respectively. This modulated equilibrium introduces resonant coupling 
between the modes with k z = ±q/2, whose frequencies are the same and can be solved 
from Eq. (5.14). Note from Eq. (5.7) that this modulation takes effect through (c a 2 /c 2 )S 
and —( oj 2 /c 2 )rdD/dr , which are functions of x according to Eq. (5.2), we thereby have 
their linear Taylor expansions about x: 
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(5.15) 


Here, uq, Lo c io , So and Dq are the corresponding values of w, cu c j, S and D for a uniform 
static magnetic field. The dependences of S and — D on the frequency difference (w — cuq) 
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and modulation (u; c j — w c *o) have been separated. By introducing variables: 
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UJQ 


Ul — UJ 0 
UJQ ' 


/a s cos(qz) = - x-^ln (w 2 S) 


(5.16) 


^ci bJciO 

UciO ' 
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we rewrite Eq. (5.15) as: 
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(5.17) 


Assuming [(w c j — w c io)/w c io] = ecos (qz) and substituting Eq. (5.17) into Eq. (5.7) result 
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Letting 'L (r) be the eigenfunction of the wave equation for the equilibrium magnetic field 
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(5.19) 


we multiply Eq. (5.18) by \H(r) and integrate over radius to eliminate the equilibrium 


terms via integration by parts, 

0 = / {yT<S 0 [/AS cos(qz) + Ps] + 
- / [/j d cos(qz) + Pd] 

+ f r Jr {“(iyS'o [/XS cos(gz) + Ps] 
We try the following form of 
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(5.21) 


with F(z) and G(z) slow functions of 2 compared with cos (qz), and k the wave number 


in the modulated field. Separating the longitudinal Fourier harmonics in Eq. (5.20), we 
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obtain two coupled equations for F and G: 
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Introducing the following integrals of the eigenfunction: 
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we rewrite Eq. (5.22) and Eq. (5.23) as: 


(5.22) 


(5.23) 


(5.24) 


F (hPs — hn + I3PD + FFs ~ hn) = — — (hsh + PDh + Psh) , (5.25) 

G (F Fs + Fk + FFd + FFs + FF = — {psh + /^dF + MsF) • (5.26) 


Recalling Eq. (5.16), we take the product of Eq. (|5.25|) and Eq. (|5.26|), and obtain 
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Multiplying Eq. (5.19) by 'l'(r) and integrating over radius lead to 
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We then rewrite Eq. (5.27) as 
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Equation (5.29) is the dispersion relation of SAW propagating in a plasma cylinder with 


periodic static magnetic field and arbitrary density profiles. 

This dispersion relation can be evaluated with a specific choice of density profile. We 


investigate a step-like radial profile of plasma density shown in Eq. (5.8) to avoid the 


continuum damping of SAW, which may occur if plasma density changes continuously in 
radius. Moreover, a realistic density profile gives multiple dispersion curves, the feature 
of coupling to other modes, which makes it difficult to study the propagating feature of 
SAW. The aim is to acquire an insight into the problem first. Recalling x = w/w c j, we 
have: 
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Through Eq. (5.21), the eigenfunction T corresponding to the step-like density profile is 
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and so we evaluate 
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Substituting Eq. (5.31) and Eq. (5.36) into Eq. (5.29), we get the dispersion relation of 


SAW for a step-like radial profile of plasma density 
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When \uo — coq\/ujo < e/2, k is imaginary and waves are evanescent along the plasma 
column, forming a spectral gap. The lower (cu_) and upper (cu+) edges of the gap can be 


calculated by setting k = 0 in Eq. (5.37): 


OJ- = COq (1 - |), co + = w 0 (l + |). 


(5.38) 


The gap width is thereby A co = — co- = eooo, which is the same to that of RLH waves 


(see Eq. 4.19). When \co — u)o\/uoo > e/2, k is real and waves are propagable along the 
plasma column, forming the spectral continuum. 


5.2.4 Defective magnetic field: gap eigenmode 

Now we introduce a defect to the system’s periodicity, similar to that employed in Chap¬ 
ter [4j to create the gap eigenmode of SAW inside the spectral gap indicated by Eq. (5.29) 


and Eq. (|5.37|). We assume that the gap eigenmode will decay away from the defect 
E v = F(z)^{r)e iqz / 2 ~ Xz + G(z)*{r)e~ iqz / 2 - Xz , (5.39) 


location as 


with A the decay constant. Similarly to Eqs. (5.21)-(5.26) but replacing k with iX, we get: 


F (FPs — i^h + FPd + hPs ~ 2 A/ 5 ) = — — (psh + PdI 3 + Psli ), (5.40) 

G (PP S + iXh + I 3 PD + hPs + 2 A/ 5 ) = — — (psh + PdH + Psh) ■ (5-41) 

The system is completed with a boundary condition for E^. An odd-parity gap eigen- 
rnode can be formed by applying the boundary condition E^(r,zo) = 0, which gives 


F = — Gexp(—iqzo). The sum and difference of Eq. (5.40) and Eq. (5.41) are: 


hPs + FPd + hPs = — 2 (hsh + PdH + Psh ) cos(qzo), 


(5.42) 
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(5.43) 


respectively. By using Eq. (5.16) and Eq. (5.28), we then solve for the frequency and decay 


constant of the odd-parity gap eigenmode: 
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Similarly, an even-parity gap eigenmode can be formed by applying the boundary condition 
E'(r, zo) = 0 (with ' over z), which gives F = Gexp(—iqzo). The corresponding frequency 
and decay constant of the even-parity gap eigenmode are: 
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(5.47) 


By using Eq. (5.31) and Eq. (5.36), Eq. (5.45) and Eq. (5.47) can be evaluated for the 
step-like radial profile of plasma density shown in Eq. (5.8): 


A = -gy sin(gz 0 ), 


A = — -g 7 sin(gz 0 ), 

for the odd-parity and even-parity gap eigenmodes, respectively. 


(5.48) 

(5.49) 


5.3 Numerical results 


To validate the theoretical analysis in Sec. 5.2, numerical computations are conducted in 
this section. We use the EMS code |124] to model the spectral gap and gap eigenmode of 
SAW. Figure[5T] shows the computational domain, which is similar to that in Chapter[4]but 
with glass tube removed. This will, for the step-like density profile employed below, ensure 
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Figure 5.1: Computational domain. The solid bar denotes a left-hand half-turn helical antenna. 
The dot-dashed line is the machine and coordinate system axis (r = 0). The coordinate system (r, 
ip, z) is right-handed with p. The double-dot-dashed line labels the inner radius of the antenna. 

that there is only one density discontinuity in the whole chamber, and it is the density 
jump constructed. The enclosing chamber is assumed to be ideally conducting, and the RF 
antenna is placed inside the plasma. To well couple SAW, which are left-hand circularly 
polarised, the RF antenna is twisted in the left-hand direction to the static magnetic field, 
and it is half-turn helical to mainly excite non-axisymmetric modes, as considered in the 
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analysis. We choose step-like radial profiles of plasma density to precisely compare with 


the analytical results in Sec. 5.2 and eliminate potential continuum damping of SAW. 


Illustrative results for broadened density profiles will be given in Sec. 5.4.1 


For a step-like-peak density profile, as shown in Fig. |5.2[ a), the discontinuity is chosen 
at ro = 0.087 m so as to ensure the total number of particles in each volume is the same. 
The on-axis plasma density is chosen to be 10 20 m~ 3 , which is sufficiently dense as assumed 



X (o)/(O c i) r ( m ) 


Figure 5.2: Plasma density profile and computed wave field in a straight plasma cylinder: (a) 
radial profile of unperturbed plasma density, (b) azimuthal magnetic field of the m = — 1 mode 
for an illustrative frequency at x = 0.64, (c) computed dispersion relation for the m = —1 mode 
(dots) and analytical result from Eq. ( 5.14| ) (dashed), together with the ratio between them (inset 
figure), (d) radial profile of the m = —1 mode at z = 14 m for x = 0.64. Solid bar shows the 
antenna location, and the endplates are located at z = 0 m and z = 30 m. 


for Eq. (5.1). We consider a single-ionised helium plasma, as did by Zhang et ah, m so 


that the on-axis electron-ion collision frequency is z/ e j(0) = 4 x 10° s . The wave field 
structure and dispersion relation of SAW are studied in a plasma cylinder with uniform 
static magnetic field. Figure [5~2j b) and Fig. 5.2 d) show the axial and radial profiles of the 
computed wave field, respectively, for an illustrative frequency at x = 0.64. A preferred 


axial coupling direction can be seen clearly from Fig. 5.2 b) due to the helical antenna 


being used. Figure |5.2|(c) shows the computed dispersion relation and its comparison 


with that from Eq. (5.14). The inset figure shows the ratio between the numerical and 


analytical results. The computed data were obtained by running the EMS for various 
frequencies, for the same conditions, and calculating the dominant axial mode number 
through a Fourier decomposition. We can see that they qualitatively agree but diverge 
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more when the wave frequency approaches the ion cyclotron frequency. This may be due 
to that the computed values of k and w, which become bigger when u —> u c i, are not 


sufficiently small as assumed in the theoretical analysis (see Eq. (5.6)). 


Next, we introduce a small periodic modulation to the uniform magnetic field, \Bq[z) — 
Bo]/Bo = 0.1 cos(40z), to form a spectral gap of SAW. The periodicity q = 40 was chosen 
to locate the gap within the frequency range of 0.5 < x < 0.9 where SAW have a long decay 
length. Running EMS for the same conditions as above, we obtain a clear spectral gap in 
both axial and radial directions, as shown in Fig. 5.3 The gap width (Ax ~ 0.067) agrees 

(a) ( b) 




Figure 5.3: Surface plots of the wave field strength as a function of driving frequency: (a) in the 
axial direction (r = 0 m), (b) in the radial direction (z = 14 m). 


well with the analytical estimate from Eq. (5.38) (Aw/u; c jo = exo = 0.069), although the 


gap location (xo = 0.63) is slightly lower than that xo (k~ = 40/2) = 0.69 from Eq. (5.14). 

Section |5.2.4| indicates that an eigenmode could be formed inside the spectral gap 
if a defect is introduced to break the system’s periodicity. Moreover, the defect should 
be located at cos(qzo) = 0 in order to form the eigenmode with frequency at the gap 
centre (see Eq. (5.44) and Eq. ( 5.46[) ) and with shortest possible width (see Eq. (5.48) and 


Eq. (5.49)). Further, the gap eigenmode can be either odd-parity or even-parity around 
the defect location, depending upon the defect configurations. Figure [5~4] and Fig. |5.5| show 
the defect configurations employed and the corresponding odd-parity and even-parity gap 
eigenmodes formed. To minimise the role of collisional dissipation and thus get a sharper 
resonant peak (see Fig. |5.4| (c) and Fig. |5.5[ c)), the collision frequency has been reduced 
to 0.0005i' e j. Figure [5~4] (a) and Fig. |5.5| (a) show that the eigenmodes are standing waves 
and their wavelengths are close to twice the system’s periodicity, which is characteristic 
of Bragg’s reflection. The decay lengths indicated by the exponential envelopes agree well 


with the analytical estimates from Eq. (5.48) and Eq. (5.49). The odd-parity and even 


parity features are indicated by the axial profiles of the azimuthal electric field around 
the defect location, shown in Fig. |5.4[b) and Fig. |5.5|(b) , respectively. Figure [5~4|(d) and 


Fig. 5.5'd) give full views of the plasma response in (x, z) space with a clear eigenmode 
peak inside the spectral gap. 

We also obtained similar gap eigenmodes for a step-like-hollow density profile shown 
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Figure 5.4: Odd-parity gap eigenmode: (a) longitudinal profiles of static magnetic field (solid 
line) and on-axis RF magnetic field for x = 0.635 (dashed line), together with the decay constant 
from Eq. (5.48) (dotted); (b) longitudinal profile of E v (on-axis) at x = 0.635, where vertical dot- 
dashed line marks the defect location and solid horizontal bar marks the antenna; (c) resonance 
in the dependence of the on-axis amplitude of the RF magnetic held on driving frequency at the 
location of the defect; (d) surface plot of the on-axis wave held strength as a function of z and x. 


in Fig. |5.6[ a). The discontinuity location is moved radially outward (tq = 0.122 m) to 
keep the total number of particles equal on each side. For a uniform static magnetic field, 
illustrative axial and radial profiles of the computed wave field are shown in Fig. |5.6| (b) 
and Fig. |5.6[ d), respectively, for x = 0.65. A preferred axial coupling direction is again 
seen clearly in Fig. |5.6[ b). The computed dispersion relation is compared with that from 
Eq. (5.14) in Fig. |5.6[ c). It can be seen that the two dispersion curves match at x = 
0.665. To form a spectral gap, we modulated the uniform static magnetic field in form of 
[Bq(z) — Bo]/B q = 0.1cos(30z). We select q = 30 to locate the centre of the gap around 


the matching point by setting k z = q/2 in Eq. (5.14), and get xo = 0.667. Figure 5.7 


shows the computed spectral gap and its comparison with the analytical estimates at 


z = 12.3 m. The analytical lower and upper edges of the gap are calculated from Eq. (5.38): 
X- = xo(l — e/2) = 0.634 and x+ = a?o(l + e/2) = 0.700. We can see that numerical and 
analytical spectral gaps agree quite well, in terms of the gap location and width. The 
spatial distribution of wave energy for frequencies inside and outside of the spectral gap 
is shown in Fig. |5.7[b), confirming that wave propagations along the plasma column are 
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Figure 5.5: Even-parity gap eigenmode: (a) longitudinal profiles of static magnetic field (solid 
line) and on-axis RF magnetic field for x = 0.649 (dashed line), together with the decay constant 
from Eq. (5.49) (dotted); (b) longitudinal profile of E v (on-axis) at x = 0.649, where vertical dot- 
dashed line marks the defect location and solid horizontal bar marks the antenna; (c) resonance in 
the dependence of the on-axis amplitude of the RF magnetic field on driving frequency near the 
location of the defect; (d) surface plot of the on-axis wave field strength as a function of z and x. 


evanescent for frequencies inside the spectral gap. The gap eigenmodes were formed by 


introducing two types of defect to the system’s periodicity, as shown in Fig. 5.8 and Fig. 5.9 


for odd-parity and even-parity respectively. We can see that they are similar to those for 


the step-like-peak density profile, shown in Fig. 5.4 and Fig. 5.5 


5.4 Experimental implementation 

In this section, we comment on the possible experimental realisation of the gap eigenmode 
of SAW in a periodic magnetic field. 


5.4.1 Radial width of density jump 

The step-like radial profile of plasma density, which has been used so far for precise 


comparison with the theoretical analysis in Sec. 5.2, is not easily achievable in experiment. 
It clearly breaks down with additional physics such as particle diffusion and cyclotron 
motion, which act to broaden the radial width of density jump. We investigate the effect 
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Figure 5.6: Plasma density profile and computed wave field in a straight plasma cylinder: (a) 
radial profile of unperturbed plasma density, (b) azimuthal magnetic field of the m = — 1 mode 
for an illustrative frequency at x = 0.65, (c) computed dispersion relation for the m = — 1 mode 
(dots) and analytical result from Eq. ( 5.14| ) (dashed), together with the ratio between them (inset 
figure), (d) radial profile of the m = — 1 mode at z = 14 m for x = 0.65. Solid bar shows the 
antenna location, and the endplates are located at z = 0 m and z = 30 m. 


of this width on the gap eigenmode formation by increasing it from A r = 0 m, as used 
for Fig. |5.4[ c), to A r = 0.01 m, A r = 0.02 m and A r = 0.03 m. Figure 5.10[ a) shows the 
density profiles broadened, and Fig. 5.10| (b) shows the corresponding resonant peaks inside 
the spectral gap. It can be seen that the peak width and location remain the same. The 
peak width was found to remain the same even when the density jump width is increased 
to Ar = 0.1 cm. Therefore, broadening the density profile shown in Fig. |5.2[ a) to that 
achievable in experiment is unlikely to affect the gap eigenmode formation. This is a 
sanity check that the gap eigenmode observed here is sufficiently resolved in the absence 
of continuum damping. Figure 5.10[ c) shows the location of this gap eigenmode relative 
to the continuum damping region. The analytical dispersion relations for uniform (same 
to the dashed line in Fig. |5.2| (c)) and periodic ([Bq(z) — Bq]/Bo = 0.1 cos(40z)) magnetic 
fields are overlaid. The location of the gap eigenmode is at the joint point of the dispersion 
relation for uniform field and the middle line of the spectral gap, because of the chosen 
defect which forces cos(qzo ) to be zero in Eq. 5.44 and Eq. 5.46 The dispersion relations for 
uniform plasma densities inside and outside of the density jump (for uniform field) are also 
shown. They label the upper and lower edges of the continuum damping region, because 
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z Cm) 


Figure 5.7: Identification of the spectral gap in simulations: (a) on-axis wave field at z = 12.3 m 
as a function of driving frequency; vertical lines show analytical predictions for the gap centre (cco) 
and the lower and upper edges of the gap (x- and x_), (b) spatial distributions of wave energy for 
antenna frequencies inside and outside the spectral gap (solid bar denotes the antenna location). 


the continuum damping only occurs when the wave frequency experiences the local Alfven 
frequency, namely to = cja = k?Bn/ a/ fj,nnurun(r). [153] We can see that the gap eigenmode 
is currently outside of the continuum damping region. We further found that deepening 
the density jump did not move the gap eigenmode location into this region. To study the 
impact of continuum damping on gap eigenmode formation, spectral gaps need be formed 
in a lower frequency region, e. g. around x = 0.2 where the dispersion relation for uniform 
field is within the continuum damping region, and a very high radial resolution across the 
broadened density jump should be employed. This will surely increase the computation 
time significantly, and we leave this work for future research. 
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Figure 5.8: Odd-parity gap eigenmode: (a) longitudinal profiles of the static magnetic field (solid 
line) and the on-axis RF magnetic field for x = 0.662 (dashed line), together with the decay 
constant from Eq. 5.48 (b) longitudinal profile of E v (on-axis) at x = 0.662, where vertical dot- 


dashed line marks the defect location and solid horizontal bar marks the antenna; (c) resonance 
in the dependence of the on-axis amplitude of the RF magnetic held on driving frequency at the 
location of the defect; (d) surface plot of the on-axis wave held strength as a function of z and x. 


5.4.2 Implementation on the LAPD 


A modification to the experiments of Zhang et al. m on the LAPD, in which they observed 
spectral gaps and continua of SAW, is a possible candidate. It is essential to ensure that 
dissipative processes in the plasma do not destroy the gap-mode resonance. [45] There exist 
mainly two types of dissipative mechanisms: continuum damping and collisional damping. 
The former needs to be avoided first of all because it will affect the mode structure globally 
and even prevent the propagation of SAW. Decreasing the radial gradient of plasma density 
to narrow the continuum damping region so that the gap eigenmode is located outside of 
the region will probably help. The latter needs to be small in order to form a sharp 
resonance peak, as pointed out in Chapter [4j The collisional damping rate of Alfven 
waves (in terms of wave amplitude) is given by | 153| 


V ct ii , 


7A = 


2 


(5.50) 
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Figure 5.9: Even-parity gap eigenmode: (a) longitudinal profiles of the static magnetic field 
(solid line) and the on-axis RF magnetic field for x = 0.679 (dashed line), together with the decay 
constant from Eq. 5.38 (b) longitudinal profile of E v (on-axis) at x = 0.679, where vertical dot- 
dashed line marks the defect location and solid horizontal bar marks the antenna; (c) resonance in 
the dependence of the on-axis amplitude of the RF magnetic field on driving frequency near the 
location of the defect; (d) surface plot of the on-axis wave field strength as a function of z and x. 


where fey and k±, fry and cry are components of the wave vector and electric conductivity in 
the directions parallel and perpendicular to the static magnetic field, respectively. Because 
of k_i <C fen and cry <C cry for the SAW, we approximate the damping rate as 


7A 



2 cry 


3 uj 2 u ei 

2 UJ c i | UJ ce 


(5.51) 


Here, the relations k\\ = lo/va and cry = cr/3 = e 2 n e /3m e i' e i (in perpendicular motions, 
the slow electrons, which have small Larrnor radii, contribute less to the conductivity than 
in parallel motions[;l2]) have been used. Equation (5.51) indicates that stronger magnetic 


field leads to larger ion cyclotron frequency and thus lower damping rate. Moreover, 
the damping rate can be further reduced by lowering the wave frequency. For the gap 
eigenmode formed here (step-like-peak density profile), the damping rate on axis is about 
8.2 x 10 3 s , which is nearly half the width of the resonance peak shown in Fig. 5.4(c) 
(around 15 x 10 3 s -1 ). The difference may be due to the fact that Eq. (5.51) does not 
consider radial boundaries but an infinite slab geometry. For the conditions employed in 








































86 


Gap eigenmode of shear Alfven waves in a periodic structure 



v.-. . 

0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 


r(m) 


( 6 ) 



Jl 40 


0 . ^ 


0.60 0.62 


0.64 0.66 0.68 

X ci) 


(c) 



Figure 5.10: Radial profile of plasma density and gap eigenmode formation: (a) density jump 
with different radial widths (solid bar labels the antenna location), (b) resonant peaks inside a 
spectral gap for the density profiles shown in (a), and (c) gap eigenmode location relative to the 
continuum damping region. 


[ TO] the damping rate is about 16.3 s 1 . Thus, the collisional dissipation on the LAPD 
can be controlled very low. Also, a bigger modulation factor of the periodic magnetic 


field broadens the spectral gap (see Eq. (5.38)), and makes it easier to resolve the gap 
eigenmode with finite damping. Detailed comparison between the parameters employed 
for the present study and those in m is given in Tab. |5.1| In this chapter, we have 


decreased the spatial periodicity of the magnetic field in order to form a larger number of 
field ripples and thereby remove conducting endplate boundary conditions, and employed 
a stronger magnetic field to weaken 7 a and thus the effect of collisional damping. The 
aim is to, for this numerical study, form a clear spectral gap and a sharp gap eigenmode, 
and verify the analytical model developed here before modelling the LAPD. 

To introduce a defect to the system’s periodicity, besides powering one magnetic coil 
with extremely strong current [TO], a conducting mesh could also be inserted radially into 
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Table 5.1: Parameters employed in this chapter and on the LAPD. 


Parameters 

This chapter 

LAPD[10] 

Length of machine 

20 m 

17.56 m 

Axial periodicity 

0.16 m 

3.63 m 

No. of periods 

127 

4 

Radius of wall 

0.15 m 

0.5 nr 

HiO 

1 x 10 20 m~ 3 

1 x 10 18 m - ' 

B 0 

1 T 

0.12 T 

T e 

4.22 eV 

6 eV 

mi 

4 amu (He) 

4 amu (He) 


the plasma. The mesh should have a radial slot to prevent the formation of azimuthal 
current, and be insulated from the enclosing chamber to avoid short circuiting the plasma. 
Because the azimuthal electric field vanishes on the mesh (E,p(r, zq) = 0), this will form 
an odd-parity gap eigenmode. To form an even-parity gap eigenmode, which requires 
dE v (r, zo)/dz = 0, a conducting ring with small radial thickness could be placed co¬ 
axially with the static magnetic field. The ring should have an axial slot and be insulated 
from the chamber, for the same reasons as for the mesh. 

Therefore, by controlling the collisional damping rate, varying the modulation factor, 
and selecting the appropriate defect, the experimental implementation of the gap eigen¬ 
mode of SAW is conceivable on the LAPD. 


5.5 Conclusion 

The gap-mode analysis of RLH waves in the whistler frequency range [45 j has been gener¬ 
alised to ion cyclotron frequencies. We applied this generalised analysis to SAW, whose 
spectral gaps and continua were recently observed on a linear plasma device (LAPD),[10j 
and utilised an electromagnetic solver 121 to model the spectral gap and gap eigenmode. 
The employed conditions, such as the ratio of wave frequency and ion cyclotron frequency, 
ion species, and geometry, are close to those on the LAPD, except the radial profile of 
plasma density which is step-like and chosen for precise comparison with the analytical 
estimates. For a uniform magnetic field, the computed dispersion relation qualitatively 
agrees with that from the generalised model. However, they diverge more when the wave 
frequency is close to the ion cyclotron frequency. This may be because the computed val¬ 
ues of k and u> are no longer sufficiently small as assumed in the theoretical analysis. For a 
periodic magnetic field, a clear spectral gap is formed in both (z, x ) and (r, x ) spaces, and 
the gap width and location qualitatively agree with the analytical estimates. Moreover, 
we found a resonant peak inside the spectral gap when the field periodicity is broken with 
a local defect. The resonant peak stands for a gap eigenmode, which is a standing wave 
and with wavelength near twice the field periodicity, a characteristic of Bragg’s reflection. 
The decay length of the gap eigenmode is consistent with the analytical prediction. We 
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further explored the effect of the radial scale of density jump on the gap eigenmode for¬ 
mation, and found that the effect is ignorable. The collisional damping rate of SAW on 
the LAPD was also calculated, for the conditions employed in m, and it was found much 
smaller than that in our simulations. Therefore, the dissipative processes on the LADP 
can be controlled at a very low level and will not destroy the gap-mode resonance. We 
conclude that experimental implementation of the gap eigenmode of SAW is conceivable 
on the LAPD, providing that appropriate modulation factor and defect are selected. 



Chapter 6 


Summary and recommendations 
for future research 


The purpose of this thesis was to explore the impact of magnetic geometry on wave modes 
in cylindrical plasmas. This has been achieved by increasing the magnetic shaping on wave 
modes in helicon plasmas. A two-fluid model was applied to describing plasma-driven 
electrostatic drift waves in a uniform field. The model provided a consistent description of 
the plasma, as well as a qualitative match to the mode structure of spontaneously-driven 
electrostatic waves. The model is amongst the first to capture the physics of flow in a 
helicon discharge. The construction of a new field-pinched helicon source at the ANU 
afforded the opportunity to explore the impact of a ramp in the magnetic field on the 
waves and plasma. In the model, antenna-driven wave fields led to plasma formation 
through the dielectric tensor and electron-ion collisions. Wave fields in agreement with 
experimental data were found with a 9.5-fold increase in the Coulomb collision frequency: 
this anomalous enhancement is consistent with other works. This analysis was amongst 
the first to describe the observed variation of wave field strength and plasma density with 
field strength. Finally, it was shown that the introduction of a periodic magnetic field to 
the system can create gaps in the continuum of whistler and shear Alfven modes, while the 
addition of a defect to the system can create gap eigenmodes. The opportunity to excite 
such modes in a simple magnetic confinement device offers the tantalising prospect to 
expand the study of gap eigenmodes of Alfven waves from fusion to basic plasma physics. 

In the remaining part of the chapter a detailed summary of each topic is given, together 
with possible extensions. 


6.1 Uniform Bq : wave instabilities in a flowing plasma 

6.1.1 Summary 

In Chapter [2j a two-fluid flowing plasma model, which was developed originally to describe 
wave oscillations in a vacuum arc centrifuge, [771 ITHl fTTFj was employed to interpret 
low frequency oscillations in a RF generated linear magnetised plasma. The following 
conclusions were obtained: 

(1) Compared to the more rapidly rotating centrifuge plasma, the drift wave insta- 
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bility, unstable at larger wavelengths, has a normalised frequency much larger than the 
normalised frequency for the centrifuge. 

(2) The difference between dispersion curves is principally due to the low rotation 
speed of RF plasmas compared to those of the centrifuge. 

(3) The measured fluctuation frequency is a factor of 3.5 lower than the predicted 
frequency, and the predicted trend of oscillation frequency with field strength for inferred 
electron temperature matches the data. 

(4) The discrepancy between the measured and predicted fluctuation frequencies may 
be attributable to the limitations of the model: assumed temperature uniformity across 
plasma column and neglected effects induced by plasma fluctuations on the externally 
applied field, and the lack of measurements in the experiment: simultaneously measuring 
plasma rotation profile and probe frequency and measuring the dependency of electron 
temperature with field strength. 

(5) The measured and predicted perturbed density profiles have a single peak in the 
radial direction, indicating the perturbed mode has n = 0, and the peak position is in the 
region of maximum equilibrium density gradient. 

(6) Qualitative agreement consolidates earlier claims that the mode is a drift mode, 
driven by the density gradient of the plasma. 

6.1.2 Recommendations for future research 

Recommendations for future research regarding this topic include: 

(1) generalising the electrostatic vacuum arc centrifuge model to include magnetic field 
oscillations to describe electromagnetic perturbations in a cylindrical flowing plasma (this 
was partly done in the appendix), 

(2) including an external drive into the generalised electromagnetic model to describe 
driven modes in a flowing plasma, 

(3) including radial non-uniformity into the ion and electron temperatures to possibly 
reduce the 3.5-factor discrepancy between computed and measured wave frequencies, 

(4) and experimentally measuring the plasma rotation profile and probe frequency 
simultaneously, and the dependence of electron temperature with held strength. 

6.2 Focused Bq : wave propagations in a pinched plasma 

6.2.1 Summary 

In Chapter [3j an electromagnetic solver based on Maxwell’s equations and a cold-plasma 
dielectric tensor was employed to describe the wave phenomena observed in a cylindrical 
non-uniform helicon discharge. Here, the non-uniformity is both radial and axial: the 
plasma density is dependent on radial and axial positions, the static magnetic held varies 
with axial position, and the electron temperature is a function of radial position. The 
following conclusions were obtained: 
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(1) With an enhancement factor of 9.5 to the electron-ion Coulomb collision frequency 
to approximate the observed attenuation, and a 12% reduction in the antenna radius to 
match the amplitude of the wave field, the wave solver produced consistent wave fields 
compared to experimental data, including the axial and radial profiles of wave amplitude 
and phase. 

(2) A local minimum in the axial profiles of wave amplitude was observed both exper¬ 
imentally and numerically, agreeing with previous studies. [HU EMI ESI [73] 

(3) Mode structure of m = 1 is consistent with the left hand half-turn helical antenna 
being used. 

(4) A possible explanation for the enhanced electron-ion collision frequency has been 
offered through ion-acoustic turbulence, which can happen if the electron drift velocity 
exceeds the speed of sound in magnetised plasmas. [53] 

(5) Density increasing in proportion to the static magnetic field has little effect on RF 
absorption, while the density level near the antenna affects the wave amplitude significantly 
at all axial locations. 

(6) The axial gradient in magnetic held increases the decay length of helicon waves in 
the target region. 

(7) The relationship between plasma density, static magnetic held, and axial wave¬ 
length is consistent with a simple theory developed previously. [45] 

(8) With increased electron-ion collision frequency, the wave amplitude is lowered and 
more focused near the antenna, which is mainly because of stronger edge heating at higher 
collision frequencies that prevent more energy transported from the antenna into the core 
plasma. 

(9) The wave amplitude prohle at the enhanced electron-ion collision frequency (by 
a factor of 9.5), which agrees with experimental data, shows consistent feature with a 
previous study that the RF energy is almost all absorbed in the near region of the antenna 
rather than in the far region. m 

(10) The antiparallel feature that wave energy is larger on the opposite side of the an¬ 
tenna relative to the direction of static magnetic field has been observed both numerically 
and experimentally, and is consistent with previous observations made in uniform field 
configurations. [Til 451 133 . [53] 

6.2.2 Recommendations for future research 

Recommendations for future research regarding this topic include: 

(1) further explaining exactly how an axially non-uniform field might affect the radially 
localised helicon mode.[35] 

(2) including different azimuthal mode numbers in the glass layer and any subsequent 
coupling to the plasma at the plasma-glass interface; this will possibly increase the com¬ 
puted wave signal which is currently weaker than measurement, 

(3) identifying independent first and second radial modes that superpose to yield a 
local minimum in wave field amplitude, 
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(4) measuring the axial profile of plasma density, and the plasma density and electro¬ 
magnetic wave field with a reversed static magnetic held, 

(5) and modelling a diverged plasma configuration which may be of interest for plasma 
thruster applications. 


6.3 Rippled B 0 : gap eigenmodes of RLH waves and SAW 

6.3.1 Summary 

In Chapter [4] and Chapter [5j spectral gaps of RLH waves and SAW were investigated in a 
plasma cylinder with periodic magnetic held, and discrete eigenmodes were formed inside 
their gaps by introducing a defect to the held periodicity. The following conclusions were 
obtained: 

(1) Longitudinal modulation of the guiding magnetic held in a plasma column creates 
a spectral gap, which prohibits wave propagation along the column when the driving 
frequency of the RF antenna is in the forbidden range. 

(2) The calculated width and location of the gap are consistent with analytical esti¬ 
mates. 

(3) A discrete eigenmode can been formed inside the spectral gap by introducing a 
local defect to the periodic structure. 

(4) Both the theoretical analysis and simulations demonstrate two types of gap eigen- 
nrode in the “imperfect” system: odd-parity mode and even-parity mode, depending the 
type of the defect employed. 

(5) The gap mode is localised around the defect and represents a standing wave rather 
than travelling wave. 

(6) The distinctive feature of the gap eigenmode is a resonant peak in the plasma 
response to the antenna current. 

(7) The gap eigenmode has two characteristic spatial scales: a short inner scale which 
is nearly twice the system’s periodicity and characteristic for Bragg’s reflection, and a 
smooth envelope which depends on the modulation amplitude and scales roughly as the 
inverse width of the spectral gap. 

(8) Experimental identifications of the gap eigenmodes of RLH waves and SAW are 
conceivable on the LAPD. 

6.3.2 Recommendations for future research 

Recommendations for future research regarding this topic include: 

(1) identifying the gap eigenmodes of RLH waves and SAW in a linear plasma device, e. 
g. LAPD, which if succeed will pave the way to further studying the interaction between 
gap eigenmode and energetic particles: important for energetic particle confinement in 
fusion plasmas, 

(2) including continuum damping into the gap-mode analysis of SAW and further 
studying its effect on the gap eigenmode formation. 
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A two-fluid electromagnetic 
flowing plasma model 


Appendix A 


A.l TOEFL model 


To describe the electromagnetic perturbations in a flowing plasma cylinder, a TwO-fluid 
Electromagnetic Flowing pLasma (TOEFL) model is developed. It is a generalisation 
of the electrostatic vacuum arc centrifuge (VAC) model developed by Hole et al.[77] to 
include magnetic field oscillations. The TOEFL model comprises: 
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with A n = eARu c i/(kBT e ) the normalisation form of the magnetic vector potential A 
and /3 ct = [fi 0 e 2 Zm(0) / (k B T e )\(Ru> ci ) 2 R 2 a convenient constant. Equation ( |A.l| ) is the 
difference between the momentum equations for ions and electrons. Equation (A.2) is the 
momentum equation for electrons in which Ampere’s law 


V x V x An = Pcte li ( Ui - u e ) 


(A.5) 


has been used. Equation (A.3) is the continuity equation for ions, and Eq. (A.4) is the 
Coulomb gauge. The electric field E and magnetic field B are normalised as: 


^ eR _ _ 5A n 

E„. = t—= —Vy — 




dr 


Bn = = V x A r 

(k B Te) 


(A.6) 

(A.7) 
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A two-fluid electromagnetic flowing plasma model 


respectively. Other terms are the same as defined in Sec. 2.3.2| This model can be easily 
benchmarked against the electrostatic VAC model by choosing a time-independent A n 
whose curl equals zero. 


A.2 Steady-state solution 


In the steady state, A n o is independent on time, the ion density is ko(r) = —r 2 , and 
the ion and electron velocities are Ui 0 = (0, fOjo, Uizo) and u e o = (0, ftl e o(r), u e2 o(r)), 


respectively. It can be seen from the axial component of Eq. (A.2) that u eS o(r ) = Uig o 


(through Eq. (A.5)). Hence, there is no equilibrium axial current in the quasi-neutral 
plasma, and consequently no azimuthal magnetic field. It is thereby reasonable to assume 
an entire axial magnetic field V x A n o = (0, 0, C?(r)).[78] which gives 


[O i0 - Ueo (r)]r = - 


e f 0(~(r) 
Pet df 


(A.8) 


through Eq. (A.5). Substituting Eq. (A.8) into Eq. (A.l) and integrating over r yield 


Czir) = Czo\ 1 - 


< 4 ( 0 ) 


R 2 [n 2 0 + 2i/;(\ t + Z)]e~ 


(A.9) 


Here, Cso = eBoR 2 uj C i/(ksT e ) is the normalised externally applied magnetic field lying 


along z. The second term inside of the square root of Eq. (A.9) represents the diamagnetic 


effect caused by the plasma rotation and thermal motion. For typical VAC parameters 
described in Tab. |2.l[ Fig. |A.l| shows the variation of this effect as a function of radius and 
ion rotation frequency. It can be seen that the diamagnetic field becomes stronger when 
it is closer to the plasma core and when the ion rotation frequency is higher. However, 
overall, it is very small compared with the externally applied magnetic field (about 3.5% 


on axis). Substituting Eq. (A.9) into Eq. (A.8) gives the equilibrium electron rotation 


(a) 

1.000 

0.995 
0.990 
0.985 
^ 0.980 
0.975 
0.970 
0.965 


W 




Figure A.l: Diamagnetic effect Cz( r )/Cso as a function of: (a) r, and (b) fho- 
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frequency 


with 


^eo(^) — "f" 


CzO 


Pet x/l - ?e- f2 


^ = C ^tf 2 [O 2 0 + 2^(A T + Z)]. 


(A.10) 


(A.ll) 


Substituting Eq. (A. 10) and Eq. (A.9) into Eq. (A.2) and integrating over r give the 
equilibrium electric potential 


Clo 


f 2 


Xo(r) =Xc~r + + O i0 Cf0 j y - f(r) + ln[l + /(r)] \ , 


(A.12) 


where Xc is an arbitrary reference potential and f(r) = \J 1 — ?e _r2 . To compare with 
the electrostatic VAC model, the plasma-induced magnetic field is now negleted, namely 
V x A n o = (0, 0, Cm); and fl e o(f) and Xo(f) are solved via a similar method. This gives: 


^eo( r ) — ^io(l + ^io) + 2t/’(A T + Z), 


Xo(r) = Xc + yy ( X + n io)r 2 + yr r 2 , 


(A.13) 
(A.14) 


and so fl e o is a constant. Figure A.2[ a) shows the comparison between Eq. (A.10) and 
Eq. (A.13), and Fig. |A.2[ b) shows the comparison between Eq. (A.12) and Eq. (A. 14). 
It can be seen that the inclusion of the plasma-induced magnetic held increases the equi- 



Figure A.2: Steady-state solutions from the TOEFL model and the electrostatic VAC model: (a) 
fleo(r), (b) xo(r). 


librium electron rotation frequency about 3% on axis, but decreases it noticeably where 
r > 1.9. The equilibrium electric potentials from the TOEFL model and the VAC model 
are the same. 
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A two-fluid electromagnetic flowing plasma model 


A.3 Perturbation treatment 


To explore normal modes of a flowing plasma system, the TOEFL model is linearised with 


perturbations of the form (see Eq. (2.7)) 

C(r, f, e , z) = Co (f) + eCi(f)e‘ (me+fc2 *-" r) , 


(A.15) 


where C is a place holder for (A nf i, A n g 1 , A n zi), In, X\, and (ipm, 0,;i, Uizi), and the 
subscripts 0 and 1 label equilibrium and perturbed parts, respectively. A matrix form of 
the perturbed TOEFL model is then obtained 


/ 


M 


V 


2^-nrl (0 


(r) 

<01 ( f ) 


<01 < 

<n< 


A n zl(f) 


= N 

ijl(^) 



AT(f) 



¥¥il(r) 

<01 (r) 


fiil(f) 

<n< ) 


^ ^izl(L) } 


(A.16) 


where 

( 


Pct r 

0 

0 

-2 

■mf (r)e r 
ipZfictr 

^Zp ct 

0 


[f(T)+rf'(r)]e r 
Pctr 
imf(r) e 

Pct r 

0 

-im5 rs £ + n iQ ipZ/3 ct r 2 + [f(r) + rf'(r)]e 
1>Zp ct r ^ 

5 rs Z + imf(r) e r 

Z/3 ct r 

0 

0 

0 


ikzf(r)e r 

Pet 

0 


izO ipZp c 

ikzfjr )^ 2 

TpZ/3 ct 

S rs 

TpZP ct r 

0 


0 

\I> 

1J7Z 

0 

0 

0 

0 


0 

0 

— X 
0 


f(r)e r 

Pet 

0 

0 

-2 

f(r)e 

^ZP ct 

6rs£ 

V 7 Z Pet 

0 

0 

0 


0 

0 

0 

0 

5 rs 

4>zp ct 

0 


and 


N = 


im[f(r) — rf (r)]e r 
Wctr 2 

(m 2 +fc|r 2 )/(r)e r 

Petr' 2 

0 

SrsZltn 2 + k 2 r 2 ') — ii/tZfi c tzor 2 — im[f(r) — r/ / (r)]e f: 

tpZI3 ot r 2 2 

im5 T . s 6+( m2 + ^|^ 2 )/(^) e ^ 

ipZ0 ct r 2 

ikz& r s 

t!>zp ct r 

0 

1 


[/M+fc§r- 2 /(r)-7-/'(f ! )] ef! 

Petr 2 

imf(r)e r 

Pf* 2 _2 

zkzf (r)e r 


fict 

imS rs £ — 'il)Z fictfljor 2 -\-[f (r) + fc|r 2 /(r) —r/ / (r)]e r 
ijjZ/3 ct r 2 2 

^rs^(l + fcf — i'lf’Z/3 c t ( mCliQ+va^f 2 — imf(r)e** 

*zp ot r 2 _ 2 

-\-i/jZfi c t Qior 2 (r)e r 

t!>ZP ct r 

0 
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r 2 

mkz /(r)e 
Petr 

0 

imf ( r)e r 

^ _2 
mkzf(r)e r 
ifjZp ct r 

_ m(S rs ^kz+iipZp c tUizO ) 

'ipZpctr 

S rs —iiJjZp c t {kz u izO~\~' UJ )r^‘-\-imrf r (r)e^ 

ipZpctf 2 

0 


f(r)f'(r)e r2 

fZPct 

imfy 

r 

ikzty 

f(r)f'(r) er2 

4>' 2 z 2 Pct 

0 

0 


—izu 


ikz 


0 


0 —izor — 20<iQr 0 ^ 

0 2f^o^ —izjf 0 


0 0 


0 —izj 


0 0 

im^ rf (r) 

ipZr -iJjZ 


ikzW 
tjjZ 

o 



0 0 


rf(r) 

ij>Z 

0 


0 

0 


0 0 

im ik~ 

0 0 / 


This matrix includes eight equations with eight unknowns, and can be solved in a similar 
way as used in Chapter [2j by reducing the number of equations through substituting 
unknowns, and shooting for the dispersion relation through matching the eigenfunction 
to certain boundary conditions. To shoot for the dispersion relation, the flute mode of 
this model, which has kz = 0 and is the starting point of the shooting procedure, needs 
to be solved first. For = 0, we have reduced the matrix to two coupled third-order 
ordinary differential equations in terms of ip^n and A n f Solving the flute mode through 
a double-shooting method is in progress. Details about the method can be found in [ T54j . 













